Reweighted techniques: definition and asymptotic convergence 



m 

O 

o 

(N 

X) 
D 



a 

o 



O 
m 
o 



X3 
O 

o 



X 
J3 



Cristian Predescu 

Department of Chemistry, Brown University, Providence, Rhode Island 02912 

(Dated: February 2, 2008) 

I define and characterize tiie reweighted methods, which are techniques used in conjunction with 
the random series implementation of the Feynman-Kag formula. I prove several convergence results 
valid for all series representations and then I specialize the results for the Levy-Ciesielski and Wiener- 
Fourier series. As opposed to the partial averaging method on which they are based, the reweighted 
techniques do not involve any modification of the physical potential. Rather, the underlying idea is to 
develop some specialized constructions of the Brownian bridge that enters the Feynman-Kag formula, 
so as to simulate the partial averaging effect. For the Levy-Ciesielski series representation, I develop 
a reweighted technique which has o(l/n^) convergence for potentials having first order Sobolev 
derivatives. It is suggested that the asymptotic convergence may reach 0(l/n^) for potentials 
having second order Sobolev derivatives. The method preserves the favorable oc log(n) scaling 
for the time necessary to compute a path at a given discretization point. For the Wiener- Fourier 
series representation, the particular reweighted method designed in the present article is shown to 
have 0{l/n^) convergence if the potential has second order Sobolev derivatives. The convergence 
constant has superior dependence with the inverse temperature as compared to the partial averaging 
method for the same series. Because the expression of the convergence constant does not actually 
involve the second order derivatives of the potential, it is conjectured that the 0{l/n'^) convergence 
extends to the potentials having first order Sobolev derivatives only. 
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I. INTRODUCTION 

Numerical recipes based on the discretization of the 
Feynman-Kag formulciii^ are the most successful ways 
to account for the quantum contributions in equilib- 
rium statistical simulations for many-body systems. The 
Feynman-Kag formula represents the density matrix and 
the partition function of a canonical quantum system as 
an infinite dimensional stochastic integral of a Brown- 
ian bridge functional i^j^ This stochastic integral is then 
approximated by a sequence of finite dimensional inte- 
grals, which in turn are evaluated by Monte Carlo tech- 
niques. The random series implementations, which 
lead to some of the most effective discretization tech- 
niques, utilize the Ito-Nisio theorem^ to construct the 
Brownian bridge in the Feynman-Kag formula. 

There are two discretization techniques of practical 
importance that can be derived from a given random 
series representation: the partial averaging2i& and the 
reweighted methodsi^ As a direct method, the partial 
averaging (PA) technique features many desirable prop- 
erties, as for instance (i) good behavior at low temper- 
ature (avoiding the so-called "classical collapse")!^ (ii) 
convergence for a wide class of potentials, including those 
that have negative coulombic singularities,^ and (iii) fa- 
vorable asymptotic convergence of up to 0{l/n^) for 
smooth enough potentialsiiS 

However, a major concern in the practical applications 
of the Feynman-Kag formula is the asymptotic rate of 
convergence as measured against the number of variables 
used to parameterize the paths n as well as the compu- 
tational time necessary to evaluate the various quantities 
involved. In this respect, the main disadvantage of the 



partial averaging method is that it involves the Gaussian 
transform of the potential. This makes the PA technique 
hard to implement if the Gaussian transform is not ana- 
lytically known. 

It has been recently recognized that there is an indi- 
rect application of the partial averaging method.=^ More 
precisely, since the properties of the partial averaging 
method are mathematically well understood, one may 
try and devise methods that (asymptotically) simulate 
the effect of partial averaging, with the hope of captur- 
ing at least some of its desirable features. As such, the 
reweighted (RW) random series techniques, which will be 
rigorously defined in the following section, are designed 
to achieve superior asymptotic convergence by simulat- 
ing the partial averaging approach. By construction, the 
reweighted techniques do not involve any modifications 
of the physical potential. Rather, they are based on some 
specialized constructions of the Brownian bridge that en- 
ters the Feynman-Kag formula. This makes the RW tech- 
niques straightforward to implement yet, for any given se- 
ries representation, there are reweighted techniques that 
achieve an asymptotic convergence better than that of 
the corresponding partial averaging method. 

As opposed to the partial averaging method, the 
reweighted techniques are not uniquely defined by the 
random series representation. In fact, we talk about a 
family of methods called reweighted techniques that are 
derived from a given random series representation. They 
share the property that their asymptotic convergence is 
characterized by the same convergence theorems. Then, 
these theorems can be used to design reweighted methods 
having improved asymptotic behavior. In this article, I 
use the aforementioned convergence theorems to design 
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two special reweighted techniques based on the Levy- 
Ciesielski and Wiener-Fourier series, respectively. The 
technique based on the Levy-Ciesielski representation 
will be shown to have o(l/n^) asymptotic convergence 
for Kato-class potentials that have first order Sobolev 
derivatives.^^ I also conjecture that the asymptotic rate 
of convergence may reach 0{\/ii?) for the potentials 
that have second order Sobolev derivatives. Quite im- 
portant, the method preserves the favorable cx log(n) 
scaling for the time necessary to compute a path at a 
given discretization point. The technique based on the 
Wiener-Fourier series is proven to have 0{\/rt') asymp- 
totic convergence for potentials that have second order 
Sobolev derivatives. However, the convergence constant 
has more favorable dependence with the inverse temper- 
ature as compared to the partial averaging convergence 
constant. Moreover, the convergence constant does not 
depend upon the second order derivatives and I specu- 
late that the cubic asymptotic convergence extends to the 
class of potentials that have first order Sobolev deriva- 
tives only. The special treatment received by the Wiener- 
Fourier series is motivated by the fact that it is the op- 
timal series representation of the Brownian bridge with 
regard to the minimization of the number of variables 
used to parameterize the paths 



In this paper, our primary interest is the development 
of the mathematical ideas lying behind the reweighted 
techniques. As such, the two Levy-Ciesielski and Wiener- 
Fourier reweighted techniques developed in the paper 
may still not be the optimal ones. Rather, they serve 
as examples demonstrating that there is a large class 
of methods which have cubic convergence for sufficiently 
smooth potentials and which are direct methods in the 
sense that they do not use "effective potentials" for their 
implementation. 



The content of the article is organized as follows. In 
the next section, I give a short review of the random 
series and the partial averaging methods with the pur- 
pose of establishing the notation. Next, I define the 
reweighted methods and justify the mathematical rele- 
vance of the definition. In the first part of Section III, 
I provide a description of the main class of potentials 
for which the reweighted method is expected to converge 
to the correct Feynman-Kag result. In the second part 
of the same section, I state some important convergence 
theorems, which are valid for all series representations. 
In Section IV, I utilize the results of Section III to de- 
sign two specific reweighted methods based on the Levy- 
Ciesielski and Wiener-Fourier series, respectively. Their 
convergence properties are reinforced by numerical stud- 
ies of a simple harmonic oscillator. I conclude the article 
with Section V, where I summarize and discuss the main 
results of the paper. 



II. DEFINITIONS OF THE MAIN RANDOM 
SERIES TECHNIQUES 

A. Random series technique and partial averaging 

In this section, I shall give a short review of the random 
series technique and of the partial averaging method with 
the purpose of establishing the notation. For a more 
complete discussion, the reader should consult Refs. (0) 
and lIlOT) . The starting point is the Feynman-Kag formula 



p{x x'-(3) ^jgg^p r ^ Cv\xr{u)+aB'^\Au\, 

(1) 

where p(x, x'; (3) is the density matrix for a monodimen- 
sional canonical system characterized by the inverse tem- 
perature (3 = 1/ {kgT) and made up of identical particles 
of mass mg moving in the potential V{x). The stochas- 
tic element that appears in Eq. (Q), {Sj], m > 0}, is a 
so-called standard Brownian bridge defined as follows: 
if {Bu, u > 0} is a standard Brownian motion starting 
at zero, then the Brownian bridge is the stochastic pro- 
cess {Bill i?i = 0, < M < 1} i.e., a Brownian motion 
conditioned on i?i = 0. In this paper, we shall reserve the 
symbol E to denote the expected value (average value) of 
a certain random variable against the underlying proba- 
bility measure of the Brownian bridge To complete 
the description of Eq J^l, we set Xr{u) — x + [x' — x)u 
(called the reference path), a = (fi,^/3/mo)^/^, and let 
Pfp{x,x'] f3) denote the density matrix for a similar free 
particle. 

The generalization of the Eq. J^l to a d-dimcnsional 
system is straightforward. The symbol now denotes 
a d-dimensional standard Brownian bridge, which is a 
vector {B^ i,B'^ 2j ■ ■ ■ d) ^i^^ t;he components being 
independent standard Brownian bridges. The symbol a 
stands for the vector (ai, a2, ■ ■ ■ , crd) with components de- 
fined by af = fi^P/mo^i, while the product aB*^ is inter- 
preted as the d-dimensional vector of components (JiB^ ^. 
Finally, x and x' are points in the configuration space 
while Xriu) — X + {x' — x)u is a line in connecting 
the points x and x' . In this paper, we shall conduct the 
proofs for monodimensional systems and only state the 
easily obtainable d-dimensional results. 

The most general series representation of the Brownian 
bridge is given by the Ito-Nisio theorem)^, the statement 
of which is reproduced below. Assume given {Afe(r)}fc>i 
a system of functions on the interval [0,1], which, to- 
gether with the constant function Ao(t) = 1, makes up 
an orthonormal basis in i^[0, 1]. Let 

Afc(t) = / Afe(u)dM. 
If is the space of infinite sequences a = (ai, a2, . . .) and 
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p[a\ - n 



(2) 



k=l 



is the probability measure on Vl such that the coordi- 
nate maps a ak are independent idcnticaUy distributed 
(i.i.d.) variables with distribution probability 

fi{akeA)^^ f e-'"/^dz, (3) 



then 



(4) 



fe=i 



is equal in distribution to a standard Brownian bridge. 
Moreover, the convergence of the above series is almost 
surely uniform on the interval [0, 1]. 

Using the Ito-Nisio representation of the Brownian 
bridge, the Feynman-Kag formula ^ takes the form 



Pfp{x,x';P) 



dP[a]exp<^ - (3 V Xr{u) 
n I Jo 

oo 



du 



fc=i 



(5) 



The independence of the coordinates Ofc , which physically 
amounts to choosing those representations in which the 
kinetic energy operator is diagonal, is the key to the use 
of the partial averaging method. Denoting by ]E„ the av- 
erage over the coefficients beyond the rank n, the partial 
averaging formula reads 



Pfp{x,x';(3) 
X exp I — /3 E„ 



= / d^(ai) . . . / d^(a„) 



V 



.(m) -f cr^ afcAfc(M) du ^. (6) 



fc=i 



To make more sense of the above formula, it is conve- 
nient to use the notation introduced in Rcf. (Q) 

n oo 

'5"(a) = ^afcAfc(u) and = ^M^), 

k=l k=n+l 

which denote the nth order partial sum and and nth or- 
der tail series, respectively. By construction, the random 
variables S^{a) and B"{a) are independent and 

oo 

i?°(a)=.^afeA,H = 5::(a) + B::(a). 
fe=i 

Moreover, the variables B^^ and i?" have a joint Gaussian 
distribution of covariances 



/c— n+l A;— n+1 

oc 

and 7„(u,t) =E„(B:b;?) - ^ Ak{u)Ak{T). (7) 



Equivalently, by using the fact that S*" (a) and B" (a) are 
independent and their sum is _B"(a), one may evaluate 
the above series to be 



k=l k=l 



k=l 



k=l 



and 



^n{u,T)^E{B"^B"^)-Y^k{u)Ak{r) 

k=l 

n 

= min(u, r) - ur - E Afc(u)Afc(r). 

fc=i 

A straightforward computation of E{B°)'^, ErS?)^, and 
E{B^B^) is presented in Appendix A of Ref. IjlOl) . 

Going back to Eq. ((HJ, one inverts the order of inte- 
gration in the exponent and computes 

E„ / duV[xriu) + aBlia)] 
Jo 
fi 

duEnV[xr{u) + (jSZ{a) + aBl{a)] (8) 



Vu,n[xr{u) + cr5"(a)]dM, 



where 

yu,n{y) 



exp 



V{y + z)dz, 



with r^j(u) defined by 



(9) 



fe=l 



(10) 

In deducing Eq. (jSJ, one uses the fact that the variable 
has a Gaussian distribution of variance E„(i3")^. To 
summarize, we define the nth order partial averaging ap- 
proximation to the diagonal density matrix by the for- 
mula 



d^(ai) . . . / dp{an) 



p'n^{x,x'-!3) 
Pfp{x,x';(3) 



X exp <^ -P Vu,n -t- crE «fcAfe(u) du L (11) 



fc=i 



In the present paper, the cosine- Fourier basis {A/c(t) = 
•\/2 cos(A;7rr); fc > 1} will receive a special treatment. 
The corresponding series representation for the Brownian 
bridge is the Wiener-Fourier construction 



B'Xa) 



k=n+l 



■E«fc 

fc=i 



sin(fc7ru) 



, < u < 1. 
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This series has superior convergence properties as demon- 
strated in Refs. ^ and (|10ll . The corresponding methods 
are denoted by the acronym WFPI. 



B. Definition of tlie reweiglited metliods 

The reweighted method was first introduced by Pre- 
descu and Doll on an "intuitive" basis as a method 
to account for the effects of the tail series in a way 
that does not involve any modification of the associ- 
ated potentiali^ As opposed to the partial averaging 
method, the reweighted method is not uniquely defined 
by the random series representation on which it is based. 
Rather, for each series representation, we talk about 
a family of methods called reweighted methods, that 
are constructed according to certain guidelines. The 
original idea was to replace the collection of random 
variables {B"{a);0 < w < 1} by another collection, 
• • • ,6„);0 < u < 1}, which is supported by an 
n-dimensional underlying probability space and is chosen 
such that: 

1. The variance at the point u of • • • , 5„), de- 
noted by r'^(u), be as close as possible to r^(u). 

2. The variables S'"(ai, ■ • • , On) and • ■ ■ , ^n) be 
independent and their sum have a joint distribution 
as close to a Brownian bridge as possible. 

It is clear that these requirements are rather lax and 
they involve a lot of experimentation in order to con- 
struct a good reweighted method. However, the recent 
availability of sharp theorems on the convergence of the 
partial averaging methodiS allows us to provide a clearer 
definition of the reweighted methods. We begin our anal- 
ysis by defining the reweighted methods and then, in the 
remainder of the section, we justify this definition. 

Definition 1 A reweighted method constructed from the 
random series J2T=i ^k^k(u) is any sequence of approx- 
imations to the density matrix of the form 



Pfp{x,x';(3) 



dfj,{ai) ... dn{aqn+p) 



X exp 



f3 / V 





qn+p 



-f CT^ afeAfe(M) (12) 



k=l 



k=n+l 



where q and p are some fixed integers and where the func- 
tions An^kiu) satisfy the relation 



qn+p oo 
fc— n+1 k—n-\-l 



(13) 



In contrast to the original definition, here the equal- 
ity r^(t) = r'^(i) is enforced exactly. As such, the 
reweighted Fourier path integral method initially utilized 
by Predescu and Doll can be no longer categorized as 
a reweighted technique. However, the reweighted tech- 
nique introduced by the same authors in Ref. (0) for the 
Levy-Ciesielski series representation remains a veritable 
reweighted method. 

To justify the relevance of Definition ^ we first have 
to introduce some additional notation. We define 

qn+p 

Kia)^ Y akKk{u) (14) 

k=n+l 



and 

ln{u,T) 



E 



qn+p 

Y A„,fc(M)A„,fe(T). (15) 

fe=n+l 



E(S:)--7„(«,u). (16) 



Then, Eq. says that 

7„(m, u) = Efi3!:'^ 
Next, we define 



Ur,{x,x',(};a)^ / V[xr{u) + (jS';^[a) + aBl{a)]Au 
Jo 

and 

Xnix,x' , P;a) = pfp{x,x';P)exp[-pUn{x,x',(3;d)]. 

This functions are analogous to the variables 
Un{x,x' , (3;d), Uoo{x,x' ,(3;d), Xn{x,x' , (3;d), and 
Xoa{x, x' , 0; d) introduced in the beginning of Section III 
of Ref. (7). To save typographical space, I shall not 
redefine the latter variables here, but I refer the reader 
to the cited paper instead. In terms of the new functions, 
the n-th order reweighted technique approximation to 
the density matrix given by Eq. (|12|) is expressed as 



prix,x';(3)^E[Xnix,x',f3;d)]. 



(17) 



We now make a crucial observation which shall even- 
tually justify Definition ^ With the new notation, the 
relation © reads 



Un{x,x' , (3;a) = E„ [Uocix, x' , (3; a)] . 
However, because of Eq. we also have 



Unix,x',l3;a) 



Er, 



U„{x,x\P; a) 



(18) 



(19) 



and this last equality implies that 

Proposition 1 The partial averaging method derived 
from any reweighted technique is identical to the partial 
averaging method derived from the original random series 
representation. 
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The reader should stop and ponder that Proposition 
is not really a theorem. Rather, the underlying idea is 
to define the reweighted method such that Proposition 
be true. The value of this idea is that almost all identi- 
ties relating x' \ (3) and p^/^{x, x'; (3) are also true of 
p^{x, x'\ ff) and (^{x^ x'] (3). For instance, an applica- 
tion of Jensen's inequality shows that 

]g ^-l3Un{x,x' ,fj-a) ^ g-/3E„a„(x,x',/5;a) _ ^-f3U„{x,x',f3;a) 

from which, by taking the total expectation and then 
multiplying with pfp{x,x'; (3), we learn that 

p^'^{x,x';P)>p^^\x,x';P). (20) 

This relation is analogous to the well-known inequality 

p{x,x';P) > p^^{x,x';/3). 

More generally, any proofs that are based solely on the 
common property of the random variables B"{a) and 
B^{a) of being Gaussian distributed variables are hkely 
to extend to the reweighted technique, too. Therefore, it 
should not be surprising that the asymptotic behavior of 
the difference 



p'^'^{x,x';P)-p'^'^{x,x';p) 



(21) 



is described by theorems similar to those that character- 
ize the difference 



p{x,x';p)- p'^'^{x,x';l3) 



(22) 



Such theorems are presented in Section III and they 
can be employed for the design of reweighted methods 
that have specific asymptotic properties. More precisely, 
let us assume that the asymptotic convergence of the 
partial averaging method is 



p{x,x']l3)- p^^{x,x';P) 



Cpa{x,x';P) 



Let us also assume that we are able to devise a reweighted 
technique for which the difference (|21|l has the asymp- 
totic behavior 



pT{x,x';f3)-pl^ix,x';(3) 



Cbw{x,x';I3) 



Then the asymptotic behavior of the reweighted method 
is 

pix, x'- 13) - pT{x, x'- (3) = [p(x, x'-H) 

~pI\x, x'- /?)] - /3) - pIH^, /?)] (23) 

Cpa(x, x'] P) - Cjiw{x, x'; f3) 



whenever CpA(a;, x'; /?) ^ {x , x' ; f3) . Actually, we 
would like the difference 

CKwix,x';(3)=CMx,x';(3)-Cnw{x,x';P) (24) 



to be as close to zero as possible. In the ideal situa- 
tion that the difference is zero, the convergence of the 
reweighted technique becomes o(l/ri*) and so, even if we 
do not know its asymptotic behavior, we know that the 
new method is faster than the corresponding partial aver- 
aging method. For the cases where a perfect cancellation 
does not happen, one may still exploit the possibility that 
the constant Cbw {x , x' ; f3) might have better properties 
than its components, as for instance more favorable de- 
pendence upon temperature or upon the derivatives of 
the potential. 

To conclude the justification of the reweighted method, 
we notice that the nth order approximation utilizes in 
fact qn + p Gaussian variables a^. As such, the asymp- 
totic convergence of the reweighted method as measured 
with respect to the actual number of coefficients used to 
parameterize the paths is 



q''Cnwix,x';P) 
{qn + pY 



(25) 



i.e., the convergence order is maintained but the con- 
vergence constant becomes g^-times larger as modulus. 
This observation also explains why we forced the maxi- 
mum number of additional parameters to scale at most 
linearly with n. Finally, Eq. H25|l indicates that, as far as 
the total computational cost of the method is concerned, 
we should try and minimize the number of additional co- 
efficients, if this does not affect the asymptotic behavior 
of the resulting method or some desirable features of the 
convergence constant {x , x' ; (3) . 



III. THE ASYMPTOTIC CONVERGENCE OF 
THE REWEIGHTED TECHNIQUE 

A. The class of potentials for which the reweighted 
techniques are expected to converge 

We begin this section by considering the class of po- 
tentials for which the reweighted techniques are expected 
to converge to the correct Feynman-Kag result. As dis- 
cussed in Refs. Q) and (IQ), the Feynman-Kag formula 
is known to represent the Green's function of the corre- 
sponding Bloch equation for a large class of potentials 
called the Kato class. A Kato-class potential is defined 
as foUowsii^ Let 



9{v) 



\V\ 

HWvW 



d=l, 
-1) d = 2, 
d>3 



and define the Kato class Kd as the set of all measurable 
functions / : M"* ^ M such that 



lim sup 



x — y\\<a 



\f{y)g{x - y)\dy = 0. 



We say that / is in A'^°'^ if lof £ Kd for all bounded 
domains D CW^. Clearly, Kd C K^. We say that V{x) 
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is a Kato-class potential or that it is Kato-decomposable 
if its negative part VL = max{0,— is in while its 
positive part V+ — max{0, V} is in K]j^^. 

Since the Kato class contains almost all physically rel- 
evant potentials, we shall restrict the analysis and the 
definition of the reweighted techniques to this class of po- 
tentials. It has been shown that if the potential has also 
finite Gaussian transform, the partial averaging method 
is convergent to the Feynman-Kag formulai^ For defi- 
nitions and further discussions, we refer the reader to 
the cited bibliography. The condition that the potentials 
have finite Gaussian transform is not necessary for the 
convergence of the reweighted techniques, which are ex- 
pected to converge even for systems that are confined to a 
certain region of space (for example, a particle in a box). 
However, as opposed to the partial averaging method, the 
non-averaged methods are known not to converge for po- 
tentials having sufficiently "bad" negative singularities, 
as for instance negative coulombic singularities liii^ii^ji^ 
Such non-averaged methods include the trapezoidal Trot- 
ter discrete path integral methods and the primitive ran- 
dom series methods. Unfortunately, they also include the 
reweighted techniques. 

To identify the problems that might appear if the po- 
tential has negative singularities, it is useful to work out 
an analogy with the proof of the convergence of the par- 
tial averaging method. Theorem 1 of Ref. (0) says that 
if the potential V{x) is a Kato-class potential that has 
finite Gaussian transform, then 



lim [/„(x,x',/3;a) = C/oo(a;,a;',/3;a) 



(26) 



P-almost everywhere (probabilists say P-almost surely 
because P is a probability measure). Of course, this im- 
mediately implies 

lim Xn{x,x' , P;a) ^ Xoo{x,x' , P;a) (27) 

n — >oo 

P-almost surely. In analogy with these equations, for the 
reweighted techniques the relations 

lim Un{x,x' , (3]a) — Uoo{x,x' , (3;a) (28) 



and 



lim Xn{x,x\(3]a) ^ Xoo{x,x' ,(3;a) 



(29) 



are expected to hold P-almost surely. However, in this 
case the requirement that V{x) have finite Gaussian 
transform is no longer necessary because the functions 
Un{x,x' , (3]a) and Xn(x,x' , l3;d) no longer involve the 
Gaussian transform of the potential. In this paper, we 
shall admit without proof that Eqs. (|28|l and (|29|l are 
true for all Kato-class potentials (functions). 
However, Eq. (|29|l does not necessarily imply 



lim / X„{x,x',(3;a)dP{a) ^ / Xoo{x,x' , P;a)dP{a). 

(30) 



To give a simple example, the sequence of functions 
fn{x) — l/(n|x|) converges to zero for any x 0, thus 
Lebesgue almost everywhere. However, 



n\x\ 



da; = cx) 7^ 



Odx. 



A similar situation happens for Eq. (|30l) whenever the 
potential has coulombic negative singularities. The left- 
hand term of Eq. H30(l is -f oo for all n > 1, yet the right- 
hand term is finite [according to Theorem 5 of Ref. (f^)]- 
In the case of the partial averaging method, it happens 
that the sequence X„{x, x' , (3; a) is dominated from above 
by the sequence E„Xoo(a;, a;', /3; a), as follows from the 
Jensen's inequality. The latter sequence has the proper- 
ties 



lim ^nXoo{x,x' ,(3;a) = Xoo{,x,x' ,(3;a) (31) 



and 



E [E„Xoo(a;, x' , /3; a)] = EXoo(a;, x' , (3\ a) ~ p{x, x'; (3) < oo. 

(32) 

In this conditions, a standard theorem from measure the- 
ory, namely the dominated convergence theorem^ii guar- 
anties that 

lim / Xr,ix,x',l3;a)dP{d) ^ / X^oix, x' , f3;d)dP{d). 
Jn Jn 

By analogy, it follows that in order to enforce Eq. (j^UJ, 
we need to restrict the potential V{x) to those potentials 
for which the sequence Xn{x, x' , f3; d) is dominated by an 
appropriate sequence of random variables Yn{x,x' ,P;d) 
satisfying the relations 

lim Yn{x,x',f3;d) ^Yoo{x,x',(3;d) 

n — >oo 

P-almost surely and 



lim 



Yn{x,x',/3;a) 



Yoo{x,x',/3;a) 



for some integrable Yoo{x,x',l3;a). This is the strategy 
we adopt in proving the following proposition. 

Proposition 2 Let V{x) be a Kato-class potential such 
that e^'^^*^^) is a K^J"^ function that has finite Gaussian 
transform for all f3 > 0. Then, 



lim p^'^{x,x';/3)^pix,x';P) 

n— *oo 

for all (x, x') ew^ xw^ and (3 > 0. If in addition 



(33) 



ZcliP) 



then we also ha 




< oo, 



lim Z^^^iP) = Z{p) < Z,i{P). 



(34) 
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Proof. Let us define the random variables 



Y„{x,x',P;a) = Pfp{x,x';(3) / exp | - /3y [^^(u) 



-<jS:{-a) + aB:{a)]} 



du 



and 



Yoo{x,x',f3;a) = pfp{x,x';f]) J exp - (3V[xr{u) 

+aBlia)]}du. 
Then Jensen's inequahty guaranties that 

Xn{x,x' , (3;a) < F„(x, x', /?; a) 

and 

Xooix,x',P;a) < Yoo{x,x' ,P\a), 

respectively. Because e"^^'^^-' is a positive K^^'^ function, 
it is Kato-decomposable and replacing V{x) with e~^^''^'^ 
in Eq. (|25|l . we learn that 

lim Ynix.x' ,l3\a) ^Yo^{x,x' ,[i\a) 

n — >oo 

P-almost surely. Moreover, using Eq. pt)|) . one readily 
computes 



Yn{x,x' ,I3\ a) 
1 

dw(27r)-'^/2 



Y^{x,x' ,l3]a) = pfp{x,x';f3) 

,„-\\zf/2~r3V[x,^{u)+c,ro(,u)z] 



< oo. 



The fact that the last value is finite is guarantied by 
Theorem 4 of Ref. ^ and the fact that e-'^^(^) is a Kato- 
class function of finite Gaussian transform [again, replace 
V{x) with in Eq. (A2) of the cited reference]. 

In these conditions, Eq. (|29|l and the dominated con- 
vergence theorem imply Eqs. (|Sn)l and (|^ . Moreover, 
we have 



and 



ZW)< 



E 



Ynix,x, (3;a) dx ^ Zci{(3) 



E 



Y^{x,x,P;a) dx = Zci{P) 



respectively. As such, if Zci{(3) is finite, Eq. 
and again the dominated convergence theorem imply 
Eq. The proof is concluded. □ 

The reader should realize that the hypothesis of Propo- 
sition |21 is likely to be satisfied for many of the physical 
systems of interest for the chemical physicist. It is trivial 
to show that if a Kato-class potential is bounded from 
below, then e~^^^^'^ has finite Gaussian transform and 
lies in K^f^ fuse the fact that e~'^^*^^' is bounded from 



above by a constant and that the constant functions are 
ifd-functions] . For such potentials. Proposition |21 guar- 
anties the pointwise convergence of the density matrices, 
as shown by Eq. Moreover, the quantum partition 

function is recovered as the limit 

Z(/3) = lim 

n — 'OO 

whenever the classical partition function is finite. 



B. Asymptotic convergence of the reweighted 
techniques 

For the remainder of this paper, we shall conduct the 
proofs for monodimensional systems and only state the 
multidimensional analogues. Also, we shall assume that 
the potential V{x) is a Kato-class potential that satisfies 
the hypothesis of Proposition[5] However, for the purpose 
of studying the asymptotic convergence, we also demand 
that the potential V{x) have finite Gaussian transform, 
so that the partial averaging method is convergent too. 
This requirement is motivated by the fact that we do not 
study directly the asymptotic behavior of the quantity 

p{x,x';(3)-p^'^{x,x';(3) 

but instead, as argued in Section II. B, we study the dif- 
ference 

[p{x,x';P) - p^^{x,x';f3)] 
-[pT'{x,x';P)-p^„\x,x';f3)]. 

The convergence of the first term in the expression above 
was extensively characterized in Ref. ifiol) . Therefore, the 
main focus in the present article is on the asymptotic 
behavior of the difference 



p'^'^{x,x';f3)-p'^\x,x';P) 



(35) 



As emphasized in Section II. B, the techniques utilized 
for proving the asymptotic convergence of the partial av- 
eraging method can be straightforwardly applied for the 
difference (|35|) as long as they are solely based on the 
common property of the random variables B"{d) and 
B^{a) of being Gaussian distributed. A little algebra 
shows that 

EnXn{x,x' , (3;a) - Xn{x,x',(3;a) = Xn{x,x' ,P;a) 

X E„ i^e-P[UAx,x' ,l3;a)-U„{x,x' ,f3;a)] - l| . 

It is clear that 

lim [Unix, x', P; a) - Unix, x' , /3; a)] = 0, 

n — 'OO 

because both terms of the above difference converge to 
the common value Uooix, x' , (3; a), as shown by Eqs. H2ti|) 
and H28|l . Then, for large n, it makes sense to expand 
the exponential in a Taylor series and retain the first 
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non-vanishing positive term in the series, which is also 
the one that controls the asymptotic convergence. In 
analogy with Eq. (26) of Ref. iflfl) . we have 



EnXn{x,x' , (3; a) - Xn{x,x',(3; a) « Xn{x, x' , P; a) 



2 



Un{x, x\ P; a) - Un{x, x' , (3\ a) 



(36) 



where we used the fact that the first term of the Taylor 
expansion cancels exactly because of the identity (|19|l . 
The mathematical significance of the symbol ~ is that 
the ratio between the left-hand and the right-hand sides 
of Eq. (|36|l converges to 1 in the limit n — > oo. As dis- 
cussed in Ref. ((lO|) . Eq. is expected to be true for 
all potentials V{x) £ (^^^^(M) which are not constant. 
The space V{x) G nQ,i^(R) is the space of all potentials 
V{x) whose square has finite Gaussian transform. For 
additional information, the reader is referred to the cited 
bibliography. Taking the total expectation E in Eq. H36|) 
and remembering that X„(2;, x', a) does not depend 
upon the coefficients beyond the rank n, we obtain 

p^'^ix, x'- P) - p^^ix, x'-p)^ |x„(x, x',p- a) 

X ijn{x,x' ,P;a) -Un{x,x' ,P;a) \. (37) 



Sharper descriptions of the asymptotic behavior of the 
difference 

pT{x.x'-P)-pI\x,x'-P) 

can be obtained for the potentials that lie in the 
classes r\o,Wa'^{^) and nQWj'^(R). As discussed in 
Ref. l(Tol) . the former class comprises the potentials 
V{x) € nQ,L^(M) for which the first order Sobolev deriva- 
tives are also nQ,L^(R) functions. For the potentials in 
the latter class, the second order Sobolev derivatives are 
nQ,L^(R) functions, too. Because their proofs were based 
solely on the common property of the random variables 
B"{a) and B^{a) of b eing Gaussian distributed. Theo- 
rems 1 and 2 of Ref. |lOr) extend to our problem in a 
straightforward fashion. Remembering that the poten- 
tials considered in the present paper are Kato decom- 
posable and satisfy the hypothesis of Proposition [3 we 
have 

Theorem 1 Assume V{x) e r\aW^'^{R). Then, 
p^'^{x,x';P)-p^^{x,x';P) 

J^'duj\rUn,r)KP,::,{u,r), (38) 



where 



",(u,t) = a^E^Xnix,x' , P;a) 



X V [xr {u) + aS'i (a) + aBZ (a)] (39) 
X V [xr (r) + aS-^ (a) + aB'; (a)] } . 



// in addition 7„(u, r) > 0, then the following stronger 
result holds 



p^'^{x,x';P)-pP^{x,x'-p) 



P' 



X du dr7„(u,T)i^f'",(u,T). 

JQ Jo 



(40) 



Theorem 2 Assume V{x) £ HaW^'^iR). Then, 



p^^(x,x';P)-p^^{x,x';P)^^ 



du dT7„(u, r) 
Jo 



where 



xV"[xr{u) + aSl{a) + aBl{a)] 
xV"[xr{T) + GS'^{a) + aB'^{a)]y 



(41) 



There is one important observation we make about the 
functions K^'2,(u,t) and Q^'^i{u,t). One may replace 
the Xn{x,x' , P]a) functions appearing in Eqs. (j39|l and 
(gH) with Xn{x,x' ,P',a) or Xoc{x,x' , P;a). It makes no 
difference as to the validity of the theorems. This obser- 
vation can be justified by using the fact that jniu, r) and 
7„(m, t)^ are positive definite integral kernels. In this re- 
spect, the reader may look up the argument utilized to 
demonstrate the equivalence between Eqs. (32) and (33) 
of Ref. Ufl). 

As discussed in Section II. B, there are an infinite num- 
ber of random series techniques associated with a given 
random series representation of the Brownian bridge. 
Some of them have better asymptotic convergence, oth- 
ers are worse. However, they all share the property that 
their asymptotic convergence is described by the theo- 
rems presented in this section. Then, we may use these 
theorems to find those reweighted techniques that have 
superior asymptotic convergence. As Theorems ^ and 
12 show, the asymptotic convergence is controlled by the 
two-point expectation 7„(w,r). As a rule of thumb, the 
better methods are those for which 7„(u,t) is as close 
as possible to 7„(u,t), the two-point expectation of the 
Brownian bridge. 



IV. APPLICATIONS TO THE DESIGN OF 
OPTIMAL REWEIGHTED TECHNIQUES 

A. A reweighted Levy-Ciesielski path integral 
construction 

In Ref. (0), the authors analyzed the asymptotic con- 
vergence of the random series technique based on the 
Levy-Ciesielski series representation of the Brownian 
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bridge. They found that the asymptotic convergence of 
the corresponding partial averaging method is 0(l/n^) 
and they also proposed a reweighted technique which 
has 0(l/n^) convergence. In this section, I shall use 
the theory developed so far to design a reweighted tech- 
nique having o{l/n?) convergence on the class of poten- 
tials V{x) € naW^''^{'R). The trick consists in finding a 
particular reweighted technique for which the constants 
Cpa{x, x'; (3) and Cy(}n{x, x'; (5) appearing in Eq. are 
equal and which uses a minimal number of additional 
variables. 

The Levy-Ciesielski representation of the Brownian 



bridge is based on the so-called Haar basis which is de- 
fined as follows. For k = 1, 2, . . . and j — 1,2,..., 2^~^, 
the Haar function fk,j is defined by 

r 2('=-i)/2, [(^-l)/2^^/2'=] 

- <^ _2('=-i)/2, t G [l/2\ {I + l)/2'=] (42) 
[ 0, elsewhere, 

where I = 2j — 1. Together with /g = 1, these func- 
tions make up a complete orthonormal basis in L^{[0, 1]). 
Their primitives 



Fk^jit) = { 2(^-i)/2[(; + l)/2'= - t], te [l/2\ {I + l)/2'=] (43) 
0, elsewhere 



are called the Schauder functions. The Schauder func- 
tions resemble some "little tents" that can be obtained 
one from the other by dilatations and translations. More 
precisely, we have 



Ffc.i(u) = 2-('=-i)/2Fi4(2'=-1u) 



for fc > 1 and 



Fk,j{u) = Fk, 



J-1 



(44) 



(45) 



for k>\ and 1 < j < 2^^^ . The scaling relations above 
have to do with the fact that the original Haar wavelet 
basis makes up a multiresolution analysis of i^([0, 1]) or- 
ganized in "layers" indexed by k. If we disregard the 
factor 2'^'^"^^/^, the Schauder functions make up a pyra- 
midal structure as shown in Fig. ^ 

Let {ak,j]k = l,2,...;j = 1,2, . . . ,2'''^} be i.i.d. 
standard normal variables and define Yo{u,a) = and 



E 



ak,jFk,jiu). 



Yk{u,a) 

Then by the Ito-Nisio theorem, 

oo 



(46) 



fc=i 



is equal in distribution to a standard Brownian bridge 
and the convergence of the right hand side random series 
is uniform almost surely. The introduction of the inter- 
mediate random variables Yk{u, a) summing all variables 
corresponding to the layer k is quite natural and we shall 
only consider the n = 2*^ — 1 subsequences of the various 
methods analyzed. In these conditions, we have 



k 

E 



Yi(u,a) 



(47) 




FIG. 1: A plot of the renormalized Schauder functions for 
the layers k — 1, 2, and 3 showing the pyramidal structure. 



and 



K(«)= E ^'("'«)- 

l=k+l 



(48) 



To define the partial averaging method, besides the sum 
(jTrjl . we need to evaluate 

OO Z 

i=k+i j=i 

quantity that enters the expression of the "effective" po- 
tential 



1 



exp 



2F2(«) 



V{x + z)dz. 
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As shown in Ref. (0), 



where 



and 



u{l - 2'=w), < u < 2-'', 
0, otherwise 



,p(u) =g„,i[w- (p- l)/2'=]. 



(49) 



(50) 



(51) 



Clearly, for a given n, the functions {gn,p{u), 1 < p < 
2^^} have disjoint support. Moreover, if we extend the 
function 

jo{u,u)^E{Bl)^ =u{l-u) 
outside the interval [0, 1] by setting it to zero, we have 

5„,i(w) = 2-S(2''u,2M- (52) 

Again in Ref. (0), the partial averaging method was 
shown to have the asymptotic behavior 

Cpa{x,x';P) 



where 



2/?3 



Cy>a{x,x';I3) 

r-l 



24mo 

e/3H||yy||2g-(i-e)/3ff 



x'^d9. (53) 



We now focus on the construction of our special 
reweighted technique. We begin by altering the shape 
of the Schauder functions for the next two layers beyond 
the layer k: 



B 



p=l 

+ak+2,2pRn,p[u) . (54) 



However, these additional 3-2'' = 3n + 3 functions are 
not chosen independently, but we assume that there are 
some nonnegative functions Co(u), Lq(u), and Ro{u) of 
support [0, 1], [0, 1/2], and [1/2, 1] respectively, such that 

CnAu) = 2-'=Co {2''u) , L„,i(«) = 2-^-Lo (2'=u) , 

i?„,i(u) = 2-'=i?o(2^u), (55) 

as well as 



Ln,p{u) = i„,i [u-{p- l)/2'=] , 
i?„,p(u) = i?„,i [u- (p- l)/2'=] . 



(56) 



Let me pause a moment and explain why this spe- 
cial construction is desirable. One of the major advan- 
tages of the Levy-Ciesielski series representation is the 
fact that in order to compute the series 5" (a) at a point 
u, one needs a number of operations proportional to 
log(n-|- 1) only. This is due to the fact that for a given k, 
the Schauder functions {Fkj{u), 1 < j < 2^^} have dis- 
joint support. By the way they were constructed, for a 
given n, the set of functions {Cn,p{u), 1 < P < 2'^} have 
disjoint support too and the same is true of the fam- 
ily {Ln^u), Rn,p{u); l<p<2'^}. Therefore, in order 

to evaluate the series S^{a) + B^{a) at a point u, one 
needs to perform a number of operations proportional to 
2 -I- log(n -I- 1). This is a desirable feature which, not 
surprisingly, is also compatible with Eqs. (|51|l and l|52|l . 
In fact, using these last equations, one easily shows that 
the definitional equation H13|l is satisfied if and only if 

7o(u, u) = u{l -u)= Coiuf + Lo{uf + Ro{uf . (57) 

In these conditions, the function 7„(M,r) defined by 
Eq. l(T3|l becomes 



i{u,t) ^ ^ \Cn,p{u)Cn,p{T) + L„Au)L„At) 

+Rn,p{u)Rn,p{T) ■ 



p=l 



I leave it for the reader to use the fact that the func- 
tions Cn,p{u), Ln,p{u), and Rn,p{u) have compact sup- 
port shrinking to zero and prove 



hm , 

™/o /o lniu,r)dudT 



5{u — r) 



in the sense of distributions. Moreover, one computes 

-1/2'= 



Jo 



7„(m, r)dudr — 2 



1/2" 



Lns{u)du 




1/2' 



-2k 




-1 2 



Co(M)dM 



Cns{u)du 
Rn,i{u)du 



1 -|2 

Lo(M)du 



1 -|2 

i?o(M)du 







Since 7„(u, r) > and 2^^ = n -I- 1, an application of the 
second part of Theorem ^ shows that 



pT{x,x'-[3)-pI^{x,x'-[3) 



Ckw{x,x';/3) 
(n + l)2 ' 
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where 



2„2 



2 

1 -|2 

Lo{u)du 







1 -|2 

Co(u)du 



1 -.2 

i?o(u)du 



-e/3ff||yy||2g-(i-e)/3H 



x')de. 



We now compare the above equation with Eq. (|5^{(l and 
see that the equahty 

holds if and only if 

r-l 



1 -|2 

Co{u)du 



Lo{u)du 



i?o(u)du 



1 

12' 



(58) 



Then, by virtue of Eq. 124() . the following result holds. 

Theorem 3 ^sswme V{x) £ naW^'^{R) is a Kato-class 
potential that satisfies the hypothesis of Proposition f^l. 
Then, a Levy-Ciesielski reweighted method specified by 
the nonnegative functions Cq{u), Lq{u), and Roiu) has 
o(l/n^) convergence if and only if Eqs. | |57| ) and 
hold true. 

However, to finish the proof of the above theorem, one 
has to construct at least on example of functions Cq{u), 
La{u), and Ro{u) that satisfy Eqs. (|^ and (|5H|l . We 
start with the approximations 

Co{u) « woFi^i{u), Lq{u) « wiF2^i{u), 
and Ro{u) « wii^2,2(u)- 

where the nonnegative coefhcients Wq and Wi are yet 
to be determined. First, we approximately enforce the 
equality (|57|l through integral values 



1 .1 
u{l — u)du — Wq / Fi i{u)^du 



+wl [ [F2,i(w)^ + F2,2(w)2] du. 
Jo 

The weights of the functions Lo{u) and Ro{u) were taken 
equal for symmetry reasons. Computing the integrals, 
one obtains the equation 

2wl +wl='i, 



which implies w = wq £ [O, \/2] and wi — \/4 — 2i 
Next, we define the function 




FIG. 2: Shapes of the functions Co(m), Lq{u), and Ro{u) 
utilized for the the construction of the reweighted technique. 



and set 



-iw) 



{u) = V4^^r^{u)F2,i{u), 



and 



It is clear that the above defined functions satisfy Eq. H57|l 
for all w G [O, ^/2\ . However, Eq. H58|) is satisfied only 
for the roots of the function 



h{w) = 



C^'"\u)du 





2 



4"'^(w)dM 



1 

12' 



(59) 



if there are any. Numerical analysis shows that the equa- 
tion h{w) = has precisely one solution on the interval 
[O, \/2] . The first few digits of the solution are 



Wr = 0.62258 . 



(60) 



The functions Co(m) — Cq^''\u), Lq{u) = Lf^''\u), and 
Ra(u) = R^''\u) are plotted in Fig. [3 From them, 
one generates the functions Cn^piu), Ln p(u), and Rn.p{u) 
that make up the additional layers in the reweighted tech- 
nique by dilatations (actually contractions) and transla- 
tions, as shown by Eqs. and (jHEl- 

I conclude this section by making a quite interesting 
observation. Using arguments similar to the ones em- 
ployed in the proof of Theorem 2 of Ref. (0), it is not 
difficult to prove that if n = 2*^ — 1, then 



p'^'^ix,x';P)^ dxi 



. / dx^ pr 



f3 



rw{u) 



u{l — u) 



, u;2Fi,i (u)2 -f (4 - 2^2) [F2s{uY + ^2,2(^^)2] 



1/2 



RW 

■ Po 



(61) 



12 




n 

FIG. 3: The constant c°^' is seen to increase to a certain 
finite value, demonstrating the 0{l/rC') convergence of the 
reweighed technique for the harmonic oscillator. 



quantity 

" " Z{I3) 

where Z{(5) is the exact partition function of the har- 
monic oscillator and Z^^^{(3) is the n-th order DPI ap- 
proximation to the density matrix. For the harmonic os- 
cillator, the short-time approximation given by Eq. I|()2|) 
can be computed analytically because it is the exponen- 
tial of a quadratic expression. Then, the Trotter com- 
position shown by Eq. (|61|l can be performed by the nu- 
merical matrix multiplication method i^?'^? The constants 
c^^^ plotted in Fig. |3 were computed in atomic units for 
an harmonic oscillator of mass toq = 1 ^md frequency 
w = 1, at the inverse temperature (3 — 10. The plot 
strongly suggests that the sequence c^^^ converges to a 
finite positive constant, demonstrating the 0(1/ rv^) con- 
vergence for our special DPI and reweighted techniques. 



where 



-[al+al+al)/2 



xexp<^-(3J V[x + {x' - x)u + aicrCoiu) (62) 
+a2aLQ{u) + a3(Ti?o(u)]du|-daida2da3. 

However, for arbitrary n, the right-hand side of Eq. 161(1 
defines a new quantity which we shall denote by 
p^^^{x, x';(3). The acronym DPI stands for discrete path 
integral methods and it is commonly used for methods 
obtained by Trotter composing some short-time approx- 
imation [in our case, the one given by Eq. H62|l ]. The 
results in the present section show that the subsequence 
n = 2*^ — 1 of p^^^ {x , x' ; (3) has o(l/n^) asymptotic con- 
vergence for sufficiently smooth potentials. However, it 
is not farfetched to assume that p^^^{x, x'; (3) has in gen- 
eral o(l/n^) convergence. If true, this would be a quite 
surprising result because the conventional wisdom is that 
one cannot obtain an asymptotic convergence better than 
0(l/n^) for DPI methods without employing short-time 
approximations that depend upon the derivatives of the 
potential. 

We can test the validity of this assumption by study- 
ing the asymptotic convergence of the DPI method (and 
implicitly, that of the reweighted method) for the simple 
example of an harmonic oscillator. Let us see that The- 
orem 121 predicts that the convergence of the reweighted 
technique is strictly better than 0{l/n^) but it does not 
say how much better. The experience with the asymp- 
totic convergence of the partial averaging suggests that if 
the potential has also second order Sobolev derivatives, 
then we might gain an order of magnitude in the asymp- 
totic convergence. As such, we expect that the asymp- 
totic convergence of the DPI method for the harmonic 
oscillator is 0{l/n^). We study the convergence of the 



B. A reweighted Wiener- Fourier path integral 
construction 

As argued in Ref. (0), the Wiener- Fourier series rep- 
resentation of the Brownian bridge is expected to be 
the optimal series representation as far as the asymp- 
totic convergence of the primitive, partial averaging, and 
reweighted methods is concerned. This is clearly demon- 
strated by the asymptotic convergence of the partial aver- 
aging method for potentials having second order Sobolev 
derivatives. As shown by Theorem 4 of Ref. ((10,). this 
convergence is 



lim [p{x,x';P)~p^,^{x,x';(3)] 



37r^mo 



p{x,x';p) [V'ixf + V'ix') 



/n21 



(63) 



Jo 



^-l3eHyti2^-l3{l-e)H 



Our experience with the Levy-Ciesielski series represen- 
tation in the preceding section suggests that it might be 
possible to devise a Wiener-Fourier reweighted method 
that has o(1/tt') asymptotic convergence. However, my 
attempts of finding a reweighted method satisfying the 
relation 

CpA(a;, x'\ 13) = Cpa{x, x'; (3) 

were unsuccessful. Though it remains an open question, 
I believe that such a method does not exist. 

In this paper, I shall present a special Wiener-Fourier 
reweighted method (RW-WFPI) for which the term in 
Eq. (|63|l that depends upon the second order derivative 
eventually disappears. Clearly, this is the best thing 
to enforce short of perfect cancellation. The resulting 
method will have more favorable dependence of the con- 
vergence constant with the inverse temperature (/3^ in- 
stead of (3'^). Moreover, the convergence constant will 
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depend only on the first order derivatives of the potential 
and I speculate that the 0(l/n^) convergence extends to 
the nQW^'^(M) class of potentials. 

The reweighted method analyzed in this section is con- 
structed as follows. First, remember that the Wiener- 
Fourier series is defined by the primitives 



. , X / 2 sinffcTTu) 
Next, let q = 4 and p = in Eq. ifT^ and define 



(64) 



(65) 



where 



satisfies the equality 

7„ (u, -u) = 7„ (m, u) V u e [0, 1] . 

Therefore, the method defined by Eqs. (jEJl, (|MI), and 
(|^ is a reweighted technique derived from the Wiener- 
Fourier series representation. While not the only one, the 
particular RW-WFPI technique that we study has supe- 
rior asymptotic convergence in a sense that will become 
clear from the expression of the convergence constant to 
be derived. 

We begin the convergence analysis by first considering 
the class of potentials nQWj'^(R), for which we antici- 
pate an asymptotic convergence of the type OiXlr?). As 
such, the scope of this section is to compute the limit 



(m, r) — sin(A:7ru) sin(/c7rT) (66) 



3n 



k=n+l 



and 



ctn = I 7„(u,M)du=— 

^0 ,,fr^ 



1 _ 1 1^1 

k=n+l k=l 



Next, we let 



(67) 



2ar, 



An,fe(w) — \ — — r„(u) sin(A;7ru), ti > 1, n <k < An. 
V oTl 

(68) 

Then, 



hm [pT{x.x'-[3)-pI\x,x';P)\ 

n — ^oo 

with the help of Theorem|21 Unfortunately, this theorem 
is not given in a particularly useful form because the 
computation of the limit 



lim / du/ dr7„(u,T)i^f';,(w,r) (69) 

is quite involved. As discussed in the preceding sec- 
tion, we can replace the function Xn{x,x' , (3;a) appear- 
ing m the definition of K^'",{u,t) [see Eq. ^] with 
Xn{x,x' , (3]d) without changing the asymptotic behav- 
ior of 



4n 



y„(u,T)= ^ A„_fc(u)A„,fc(r) = ^^r„(M)r„(r) 



k=n+l 



3n 



4n 



X y^ sin(fc7r?i) sin(fc7rr) 

k=n+l 



1 y.1 

du dr7„(w, r)ii:f'",(u,T). 
Jq 



Doing so and noticing that [see Appendix C of Rcf. (lO-) ] 



E {X„(x, x', P; a)En [v'[xr{u) + aS^^ia) + aB'^id)]V'[xr{T) + aS'^id) + aB'^id)]}} 

oo 2k 



we obtain 

/ du dr7„(u,r)<-(u,r)«a2/ / ^ 

JQ Jo JO JO j.^0 

X7„(u, T)'=+iE{x„(a;, x' , P; d)vl'^^^ [xr{u) + aS^id)] 
xv'-^^t'\xr{T)+aS'^{d)]'jdudT. (71) 
Because the kernels 7„(u,t)''' are positive definite for all 



A; > 1, the terms of the series appearing in Eq. ifTHl are 
positive. We now argue that we can truncate the series 
to the first two terms without changing the asymptotic 
behavior, at least as far as the 0{l/n^) convergence is 
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concerned. That is, we want to prove that We have the estimate 

"fcT 



"'0 Jo Jo Jo j.^Q 

X7„(u, t)'=+1e{x„(x, a;', /3; a)^!';^'^ [xr{u) + aS^ia)] 
xFi';+''[x,(r) + aS'^ia)]'jdudT + o (l/n^) . (72) 
I 

OO 2k rl rl 

^^'E^/ duy^ dT7„(^,r)'=+2E{x„(a;,a:^/?;a)Fi'+''K(^.) + a^:(a)]Fi'+^^ (73) 



o-^ / du dT7„(w, r)^S„(u, r), 
Jo Jo 



where the function 



S„(u,t) = E 



X. 



n{x,x',p; a) (^E„ {i^"[x,(m) + aS^ia) + aB';:{a)]V"[xr{T) + <jS!^{a) + <yB!^{a)]} 
{E,X'[xr{u) + a^:(a) + ^K(a)]} {E„r'[x,(T) + aS:^{a) + aK(a)]} 



has the property 



Hm S„(m, t) = 0. 



Using this relation together with Eq. (|B1I) . we deduce 
that the last term of Eq. (|7S|) has the property 

lim 71^(7^ du dr7„(u, t)^S„(u, r) 
Jo Jo 



— ^ hm / ^„(u, u)du = 

cJtt^ n— ►oo 



and therefore, Eq. (|72|) is proved. 

From Eq. (|73 and again Eq. IBip . we learn that 

lim du dr7„(u, T)i^f '",(m, r) = /(x, x'; ^5) 

Jo Jo 

+^ E{Xoo(a:,x',/3;a)F"[a:,(u) + aB0(a)]2}dM 



(74) 



= I{x,x';(3) + J p{x,y;u/3) 

xp[y,x';il-u)p]V"{yfdu, 



where 



I{x, x'; (3) — lim 



1 /.I 



7ti(?^,T") 



/o Jo 

X e| X„ (x, x' , /3; a)Fl_„ [x^ (u) + cr5;' (a) 
xK ,.[a;^(T) + f7S';'(a)]}dudT. 



Now, let us consider the function 



-^^.iu, t) - a2E{x^(x, x', a)T/'[x,(u) + aBO(a)] 

xF'[x^(t) +f7B^(a)]}. (75) 



[this is Eq. (35) ofRef. fW)]. Proceeding in a way similar 
to the case of the related functions K^'^,{u, r), one may 
replace Xr^{x,x' , (3;d) with Xn{x,x' ,/3;a), use the fact 
that 7n(u, T)7n{'u, t)*' are positive definite kernels for all 
fc > 0, and employ Eq. ljB2|l of Appendix B to justify the 
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equality 



1 r-l 

Jo JQ Jo 



1 /.I 1 „2k 
2 1 / ^ 

^ kl 

fe=0 



X7n(w, t)7„(", T)^-E{x„(a:, a;', /3; a)!/!';;^'^ iM^) (76) 
+t75^(a)]Fi';^'^[:r,(T) + a5;?(a)]}d?/dT + o (l/n^) . 
From Eqs. (|76|l and Hij2|l . one obtains 



\im du dT^n{u,T)K^ ,j.,{u,t) ^ I{x,x' ; (3) 



r-l rl 

a;, a:' 

'"^ Jo Jo 

+ P(a;,y;"/3)p[y,a:';(l-u)/3]y"(y)2. (77) 

We eliminate the unknown fimctions I{x,x'][3) from 
Eqs. Gill and lO to obtain 

lim / du I dr7„(u, T)i<Cf''",(w, r) 
Jo Jo 

= lim rv^ I du( dr7„(u,T)iff^,(w,r) (78) 
Jo Jo 

4 



cr 



y; y-P)p[y, x] (i - u)/3]y"(y)^ 



With the help of the this last equation, Theorem [3 and 
Eq. HBll) . we learn that 

lim n3 [pT{x,x'-p)~pl^{x,x'-f3)\ 

lim [ duf dr7„(w,T)xf^^,(u,T) (79) 



2 n- 
2^4 /•! 



Jo 



p{x,v-up)p[y,x'-{l-u)l3]V"{yf 







The remainder of this section deals with the limit 

r-l fl 



lim du dT^n{u,T)K'^^,{u,T). (80) 

Jo Jo 

As discussed in Section IV. B of Ref. l(lol) . the integral 
appearing in this limit can be written as 

1 j*U 

du dTJn{u,T)K{u,T), 



where 



Jo 



K{u,T)^K!^^^,{u,T)+Kf^,Ju,T). 
The limit in Eq. H8UI) becomes 



lim 

n — >oo 3 



4n „1 „M 

/ dw / dTii:(M, r)r„(u)r„(T) 
, 1 Jo Jo 



fe=n+l 



271 /'"'^ /"^ 

X sin(fc7ru) sin(fc7rr) = lim \J / du d 



/£=ri+l ' 
xK{u,T)fn,k{u)fn,k{T), 



where we used the fact that the product nan converges 
to I/tt^ [see Eq. (|A9|I of Appendix A] and introduced the 
functions 

fn,k{^) = rn[u) sin(fc7ru) 

[see Eq. (|A1|) of Appendix A] . 

In these conditions, Theorem of Appendix A shows 
that the hmit in Eq. l|Sn|) is 



-l^\j^\l + H{a)fda\[K{m + K 



(1,1)] 



1 d 



-^[i.(l,l)-WO)] + ^y^ -Kiu^r) 



du, 



where the function H{a) is defined by Eq. ijAMJ. Now, 
Eq. (68) of Ref. Hfl) says 

K{0, 0) = K{1, 1) = a^p{x, x'- (3) [V'{x)' + r{x')^] , 

while Eq. (72) of the same reference shows that 



' d_ 





K{u,t] 



du 



mo 



X / / Pix,y;ul3)p[y,x';{l~u)P]V"{yfdydu. 



The limit in Eq. H8U|) becomes 

^1 [1 + H{a)fda\pix,x';p) [V {xf + V {x'f 



1/4 



+ 



/ p{x,y;u/3)p[y,x';{l-u)p]V"{yfdydu. 
47r4mo Jo Jr 

Replacing the last expression in Eq. H79|l and perform- 
ing some straightforward simplifications, we conclude 



lim [p^'^ix,x';P)-p^„\x,x';p)] 



37r*mo 



X I y ^ [1 + Hia)fda \ p{x, x' ■ (3) [V'{xf + V'{x'f] 



+ 



mrhnl Jo 



-l3eHyii2^-l3{l-e)H 



')de. 



Combining this result with Eq. (|63|l . we end up with the 
following theorem. 

Theorem 4 Assume V{x) G naW^'2(M) is a Kato- 
class potentials satisfying the hypothesis of Proposition\^ 
Then, 



lim n' [pix,x';P)-p^'^ix,x';P)] = 



NoTT'^mo 

xpix,x';p)[V'ixy + V'{xr], (81) 



where 



No = 3 - /^Jl + ff(a)]'da| 



(82) 
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The integral in Eq. (|82() has been computed numerically, 
the resulting value of the constant Nq being 



iVo ~ 15.0045. 



(83) 



In the case of a d-dimensional system, the statement 
of Theorem 0] becomes 



. i=l 



[d,v{x)f + [d,v{x')f 



(84) 



where diV{x) denotes the first order partial derivative 
dV{x)/dxi. For a better understanding of this formula, 
the reader is advised to compare it with Eq. (76) of 
Ref. 113). 

We conclude this section by numerically verifying the 
findings of Theorem^lfor the simple case of an harmonic 
oscillator. In fact, we shall verify the following corollary, 
which is a direct consequence of Theorem 0] 



Corollary 1 



lim n 

n — >oo 



; Z{(3) - Z^l^^iP) 4 pix; f3)AZ{x; /3)dx 



Im pix; f3)dx 



whe 



AZ{x;P) = ^^rixf. 

Corollary 1 gives an estimate for the relative error of the 
partition function as an average of a suitable estimating 
function. In practice, such an average can be evaluated 
by Monte Carlo integration if so desired. 

The harmonic oscillator we study is described by the 
quadratic potential V(x) = mi^tjp'x^ jl. The computa- 
tions are performed in atomic units for a particle of mass 
Too = 1 and for the frequency lo = 1. The inverse tem- 
perature is set to /3 = 10. Since the density matrix for 
an harmonic oscillator is analytically known the con- 
vergence constant for the relative error of the partition 
function can be evaluated directly from Corollary 1 to be 



crw 



f^p{x;f3)AZ{x;(3)dx 
/r pix; /9)dx 



= 0.684. 



I remind the reader that, according to Eq. (|25|l . the actual 
convergence constant as measured against the number of 
variables used to parameterize the paths is 4"^ • 0.684 w 
43.8. 

We can also evaluate the convergence constant crw by 
studying the limit of the sequence 



„RW 



, zif3)-zr{P) 



The computation of the terms Z^ {(3) is discussed in 
Appendix C. Because the even and the odd sequences 




FIG. 4: Both the even sequence C2n and the odd sequence 
C2n+i are seen to converge to crw = 0.684, which is the value 
predicted by Corollary 1. 



and 



have a slightly different asymptotic be- 



havior that generates an oscillatory pattern in plots, we 
plot them separately. In these conditions. Corollary 1 
says that 



1- RW 

Crw = hm c„ 

n — >QO 

and the prediction is well verified for the harmonic oscil- 
lator, as Fig 21 shows. I regard this numerical example 
as strong evidence that the analysis performed in the 
present paper is correct. It is hard to believe that the 
exact agreement between theory and numerical analysis 
for the harmonic oscillator is accidental, especially in the 
view of the complicated definitional expression for the 
constant iVo [see Eq. 



V. SUMMARY AND DISCUSSION 

In this paper, I have defined and characterized the 
reweighted methods, which are computational techniques 
used in conjunction with the random series implementa- 
tion of path integral simulations. I have also analyzed the 
class of potentials for which the methods are expected to 
converge and I stated some general convergence theorems 
of use for arbitrary basis. Using the developed theory, I 
have constructed two special reweighted techniques start- 
ing with the Levy-Ciesielski and Wiener-Fourier basis, 
respectively. The methods were shown to have superior 
asymptotic behavior as follows. The reweighted Levy- 
Ciesielski method has o(l/n^) asymptotic convergence 
for potentials having first order Sobolcv derivatives and 
it preserves the favorable oc log(7i) scaling for the time 
necessary to evaluate a path at a given discretization 
point. For potentials having second order Sobolev deriva- 
tives, I conjecture that the asymptotic convergence of 
the reweighted Levy-Ciesielski method is 0{l/n^). The 
reweighted Wiener-Fourier method was shown to have 
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0(l/n^) convergence for potentials having second order 
Sobolev derivatives. The convergence constant has more 
favorable dependence with the temperature and with the 
derivatives of the potentials as compared to the constant 
for the corresponding partial averaging method. Since 
the expression of the convergence constant does not ac- 
tually depend on the second order derivatives of the po- 
tentials, it is suggested that the 0{l/n^) convergence ex- 
tends to the class of potentials having first order deriva- 
tives only. 

The examples analyzed in the present paper sug- 
gest that in general for any basis one may construct a 
reweighted technique which has an asymptotic conver- 
gence equal to or faster than that of the partial averaging 
method. More important, this fast asymptotic conver- 
gence is achieved with the help of some special construc- 
tions of the Brownian bridge that enters the Feynman- 
Kag formula and not by modifying the potentials. This 
make the numerical algorithms fast and straightforward 
to implement. The main conclusion of this paper is that 
there are path integral techniques capable of reaching 
0{l/n^) asymptotic convergence without any recourse 
to effective potentials. 
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APPENDIX A 

In Section IV, we introduced the functions 
fn,k(u) = rn{u) sm{kTTu), n>l, n < k < 4n, (Al) 
where 



rn(u) 



/7„(u, u) 



(A2) 



The functions 7„(u, r) and 7^(u, r) are defined as follows: 

2 ■^S^ sa\{k'Ku) sinf/cTrr) 
7„(u,t)-— }^ — 



k=n+l 

min(?i, t) ~ UT 



fc2 



2 ■^-^ sin(A;7ru) sin(A;7rT) 

2 (^3) 



k=l 



and 



°(u,t) = -— ^ sin(A:7rM) sin(A:7rT), (A4) 



4n 



3n 



respectively. The sequence a„ appearing in the previous 
formula is 

•"J k=n+l fc=l 

(A5) 

The purpose of this appendix is to study the asymp- 
totic properties of the functions fn,k{u), which are used 
in the construction of the reweighted method. In par- 
ticular, we are interested in establishing an integration 
by parts formula. Before giving the main results of the 
section, it is convenient to rearrange the expression for 
the functions r„(u) as follows. We write 



2 sin{kT:u)^ 1 



oc ^ 
E l2 



k=n+l 



1 cos(2fc7ru) 

k=n+l 



where 



9n{u) = — E 

Similarly, we have 



fc=n+l 

= Q!„[l - (A6) 
cos(2fc7ru) 



k—n+l 



2a, 
3n 



l°n{u,u) = ^ ^ sin(/c7ru)2 = a„[l (AT) 



k=n+l 



where 



^ in 

°n{u) = ^ E cos(2fc7rM). 



Therefore, 



1 - gnju) 



(A8) 



Using Eq. ljA8|l . one can prove the following proposi- 
tion. 

Proposition 3 rn{u) — > 1 on (0,1). 

Proof. By letting n — )■ oo in the sequence of inequalities 

_ /"OO ^ -1 

n / -2 1 1 
— n X ax < n y —77 

^"+1 " 



< n X dx — 1, 

In+l 



one readily proves 



k=n+l 



lim n7r^a„ = 1 and -r- < a„ < -^r—- (A9) 

n^oo 7r^(ri+l) TT^n 



18 



Now, 



cos(2fc7ru) 1 
sm{TTu) ^ = - 2^ 



k=n+l 



sin[(2fc + 1)7 

fc2 



1 ^ sin[(2fc - l)7ru] _ 1 sin[(2?i + 1)™] 



fe=n+l 



k=n+l 

and therefore, 

00 

sin(7ru) 

00 

+5 E 



fc2 (fc+l)2 



(n + 1)2 
sin[(2fc + l)7ru] 



cos(2fc7ru) 



k=n+l 



1 1 
< -- 



1 



1 



P (fc+l)2 



2 (n + 1)2 

1 



From Eqs. (|A9|I and IjAlOjl we learn that 
1 1 



(n + l)2 



\9n{u)\ < 



1 1 
< -- 



(AlO) 



(All) 



7r2Q;„(n + 1)2 sin(7ru) n sin(7rM) 
On the other hand, we have 

sin(7rw) cos(2fc7rM) 

k=n+l 
^ in 

= - ^ I sm[(2A; + 1)ttu\ - s\n[{2k - 1)ttu\ | 



k=n+l 



= - {sin[(8n + 1)ttu] - sm[(2?i + 1)ttu]} 



and therefore. 



in 



cos{2k-Ku) 



k=n+l 

The last inequality implies 

\9°n{u)\ < 



< 



1 



sin(7rM) 



1 1 

3n sin(7ru) 



(A12) 



for all u e (0, 1). Equations I|A11|I and l|A12|l readily 
imply 



lim r„(u) = 1, VmG (0,1) 



(A13) 



and the proof of the proposition is concluded. □ 
Observation. The technique utilized in establishing the 
bounds given by Eqs. (|ITl1l and l|IT2)l can be used to 
establish a similar bound for the series 



E 

k=n+l 



sin(2fc7ru) 



which will be needed later. More precisely, we have 
sin(2/c7ru) 1 cos[(2fc — 1)ttu] 



sm 



(ttu) 



k=n+l 



k=n+l 



1 cos[(2fc + 1)ttu] 1 cos[(2n + 1)t:u] 

9 2^ 



2 ^ k 

k=n+l 

00 - 

2 ^ 

fc=n+l 

This equality implies 



fc+ 1 



n + 1 



cos[(2fc + l)7ru] 



sin(7ru) 



E 

fc=n+l 

00 



sin(2fc7ru) 



1 1 
< -- 



fc=ri+l 



/C+ 1 



2n+ 1 
1 



?i + 1 ' 



from which 



E 



sin(2A:7ru) 



k 



< 



n + 1 sin(7ru) 



(A14) 



The convergence of r„(u) to 1 is not uniform. In fact, 
the functions r„(u) present singularities of the form u~^^^ 
at the end points and 1. This can be proved as follows. 
From Eq. HA3|I we obtain for u small 



7„(m, u) = u-u'^ — - 



sin(A:7ru)^ 



u. 



k=l 



On the other hand, from Eq. (|A4I) . we have 



lniu,u) 



2an 
3n 



sin(fc 



nu) 



k=n+l 

and from Eq. HA2|I . we deduce 

/ in 



2a 
3n 



n 2 
TT 



E 

\k=n+l 



rn{u) 



2a,y ^ fc' /(37i) 



^fc=n+l 



'1/2 



,-1/2 



for u small enough. The situation is similar at the other 
end point and follows by symmetry [trivially, r'„(u) = 
r„(l — u)]. However, the functions fn,k{u) themselves 
have no discontinuities because they clearly behave like 
1*1/2 jjgg^j. ^j^g gj^j^ point and like (1 — m)1/2 near the end 
point 1. 

The next few propositions in this appendix help es- 
tablish Proposition|Slwhich characterizes the asymptotic 
behavior of the primitives 



Gn,k{t) = / [r„(u) - 1] sin(A:7ru)du. 
Jo 

It is also helpful to introduce the functions 

(l>n,k{u) = [gn{u) - 5°(u)] sin(/c7ru) 



(A15) 
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and 



ipn,k{u) = |r'n(M) - 1 + 2 [9niu) - g^iu)] j siii{k7Tu). 
The primitives of these two functions will be denoted by 

^n,k{t) = / (f>n,k{u)du and ^'„,fe(i) = / V«.fc(u)du, 
Jo Jo 

respectively. Therefore, we have 

Gn,k{t) = -\'^n.k{t) + fn.fcW. (Al6) 

Finally, let us notice that 



fn,k{u) = sin{k'Ku) - -(j)n,k{u) + i^n,k{u) (Al7) 

and that the primitives of sin(fc7ru) are known analyti- 
cally. 

Symmetry arguments which shall be invoked at the 
end of the appendix show that it is enough to establish 
the asymptotic behavior of these primitives on the in- 
terval (0, 1/2] and therefore, we shall turn our attention 
to this problem. We begin by analyzing the behavior of 
the primitives ^n,k{t) on the interval (0, 1/2] and for this 
purpose, we first need to prove the following proposition. 



Proposition 4 Let qi — > oo and pi oo be two se- 
quences of natural numbers such that qi/pi — > a. Assume 
a^l and let t G (0, 1/2]. Then 



sm{qiTTu) siii{pt'n-u) 1 
lim / — — r au — — in 



1 + a 



1 -a 



We also have 



(A18) 



sup 

te(04/2] 



sm{qT:u) sm{pT:u) 
sin(7rM) 

i + q/p 



du 



<l^ln 

- 2 TT 



1 - q/p 



(A19) 



sup 

*e(o,i/2] 



sin[(2g + 1)^ 
sin(7rw) 



-du 



< 



l + -ln(2g+l^ 



(A20) 



Proof. First, let us notice that by the symmetry of the 
function 



/' 

Jo 



smlqiTTu) sm{piTru) 
sin(7ru) 



du 



in the variables qi and Pi, we may assume without loss 
of generality that a < 1. Moreover, since Eq. (|A19|) is 



trivial for qi = Pi, we may also assume that qi < pi. In 
these conditions, Pi — qi oo as i — *■ oo. Then 



sin(gi7ru) sm{pinu) 



du 



sin(7ru) 

* cos[{pi - qi)7ru\ - cos[{pi + qi)nu\ 







sin(7rM) 



du 



Jl-qi/pi 



PiSm{piTTUv)f{u)dvdu^ (A21) 



where 



/(") = 



sin(7ru) 



The function f{u) is well defined and in fact it has con- 
tinuous derivatives of any order on the set (0, 1/2]. The 
function f{u) and its first order derivative can be shown 
to be positive and increasing on the interval (0, 1/2]. Us- 
ing this observation, the reader may prove that 

sup = 7r/2 and sup = tt. (A22) 

■uG(0,l/2] tie(04/2] 

Inverting the order of integration in Eq. (jA21|l and then 
integrating by parts against the variable u, we obtain 



l+qi/pi /•* 



2 Jl-qi/pi 



Pi sm{piTTuv) f {u)dudv 



1 







1+qJPi 



1 



+ 



l-q./p. 

2 f^+qilvi ft I 
nv 



1 - cos{piTrtv)f(t)]dv (A23) 
cos{piTTuv) f' {u)dudv . 



2 Jl~qi/pi Jo 

From Eqs. (|X2Tll . (|X22|) . and (|X23|) we learn that 



sm(qiTru) sm{pi'Ku) 



< 



sm{T:u) 
i+qi/pi I 



du-Un('-±^ 
2ti \l-qi/Pi 



l-qi/pi 



■ cos(piTTtv)dv 



(A24) 



l+q./pi 



■ cos{piiTuv)dv 



du. 



Integrating by parts the first term of the right-hand side 
of Eq. HA24P , we obtain 



1+9. /p. 



i-qi/pi 



■ cos{piTTtv)dv 



■l+qi/pi Y 

< I —dv ~ In 

'i-g./p. ^ 



1 + qj/Pz 

(A25) 



and so 



sm{q.iTTu) sm{piTTu) 1 / 1 + qi/pi 

-du — — m ' 







sin(7rM) 



<i'- + lt)ln(l±^^ 

- 4 2 J \l~q^/p^ 



2ti \l-qi/Pt 

<ilnfi±^ 

2 \l-qi/pi 



. (A26) 



20 



Then, Eq. HA19|) follows trivially from Eq. (|A26|I . 

On the other hand, by integration by parts, we obtain 



1+9. /p. I 

— cos{piTTtv)dv 



3in[7r(pj + qi)t] 



cl+qi/pi 2 

, , -^sin{piT:tv)dv 



sin[7r(pi - qi)t\ 1 



< 



1 / 1 



1 



-Kt \pi + qi Pi - q, 
1 1 



1 



TTtpi 

2 



(A27) 



Thus, 



l-q,/pi l + q,/p,J TTt{p.,-qi)' 



lim 



r^+qi/Pi I 

/ — cos{pi7Ttv)dv 

Jl-qi/p, V 



= 



(A28) 



and by the dominated convergence theorem [use also 
Eq. iX25ll ]. 



lim 

i — *oo 



l+q^/pi I 

— C08{piTTUv)dv 
l-q./p, « 



du = 0. (A29) 



Eq. (|A18(I then follows easily from Eqs. l|A24p . I|A28(I . 
and HA29|1 by letting i ^ oo. 

To prove Eq. (|A20p . we first notice that 



2 sin(7rM) cos(2fc7ru) = sin[(2A: + l)7ru] 
fc=i fe=i 

q 

— ^ sin[(2/c — l)7ru] = sin[(2g + l)7ru] — sin(7ru) 



fc=i 



and so, 

sin[(2(7 + l)7rul _ „ , 
— ^ — ^ = l + 2> cos 2A:7ru . 

^ ' k=l 

Then, by the positivity of the integrand, 
sin[(2g + l)7ru]2 



sup 

te(oa/2] 



sin(7rM) 

sin[(2g+ l)7ru]^ 
Q sin(7ru) 
q ^1/2 



-du 



1/2 



du = sin[(2(7 + l)7ru]di 



+2 2^ / cos(2fc7ru) sin[(2g + l)7rM]du. 

Using the trigonometric identity 

2 sin(a) cos(/3) = sin(a — (3) + sm{a + (3) 
the above quantity becomes 

r.l/2 q 1.1/2 



I sin[(2(7 + l)7r'u]du + / sin[(2(7 - 2fc + l)7rM]du 

+ y I sin[(2g + 2A: + l)7ru]du = 

ttVo ^2q + l 

"^TT ^ 2(7 - 2/s + 1 ^ TT ^ 2(7 + 2/;;+ 1 tt ^ 2fc + 1 ' 
fc=l ^ /c=l ^ fc=0 



We write the last sum as 



and observe that 



fc+i 



-da; > 

X 



k+i 



xdx 



1 



fc + 1/2' 



by the convexity of the function \/x and by Jensen's 
inequality. Therefore, 



sup 

te(o,i/2] 



sin[(2(/ + l)7rii] 
Q sin(7ru) 



■du 



l + -ln(2(? + i; 



The proof of Eq. I|A20|1 and of the proposition is con- 
cluded. □ 



Proposition 5 Let ki ^ oo and ^ oo 6e two se- 
quences such that Ui/ki — > a. Assume a ^ 1/2 and 
t e (0,1/2]. Then 



lim [ki-K^n.MW] = 1 - aln 



1 + 2a 



1 - 2a 



1 


K 


1 + 2a 


)-"( 


1 + 8a 








)] 






1 - 2a 




1 - 8a 





(A30) 



Moreover, there are positive constants c'l, C2, C3, and C4 
such that 



sup \k'!T^n,kit)\ < c[ — C2 In 

te(o,i/2] 
if k ^ 2n + I and 



2n+ 1 



(A31) 



sup |fc7r$„,fc(i)| < cjj + C4 In {2n + 1) , (A32) 



te(oa/2] 
ifk = 2n+l. 

Proof. We begin by analyzing the term 



gn{u) sin(fc7ru)du. 



Integrating by parts twice, we obtain 
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/cTT / (7„(u) sin(fc7ru)dw : 



cos{kTTt) 
k-K 

oo 

E 

Lj=n+1 



E 

j=n+l 



cos(2j7ri) 



■ E 

2 



* cos(2j7ru) . 
^ sin(fc7rujau 



—{-{ f -\ 



OO 



■ E 

i=n+l 



sin(2j7ru) cosf/cTrt) 

cos(K7ru)au = 1 r 

i I" 



(A33) 



cos(2j7ru) 



2 sin(fc7ri) 



E 



sin(2j7rt) 



The integration term by term performed in the above is uniformly convergent. The last term in Eq. (|A33(I can 
equation is justified by the fact that the series be evaluated as follows 

cos(2j7ru) 



/ sin(2j7ru) cos(fc7ru)dM = lim / sin(2j7ru) cos(fc7ru)d 



j=n+l 



2 Jo 



sin[(2n + 1)ttu\ sin(fc7ru) 
sin(7ru) 



j=n+l 



du -\ — lim 

2 W^cxj 7Q 



sin[(2A^ + l)7ru] sin(A:7ru) 
sin(7rit) 



du. 



(A34) 



Noticing that sin(fc7ru)/sin(7ru) is integrablc on [0,1/2] 
and employing the Ricmann-Lebesgue lemma, we learn 
that 



lim 

N-^oo 



3in[(2iV + l)7ru] sin(fc7ru) 
sin(7rM) 



du = 0. 



Therefore, 



oo „t 

/ sin(2j7ru) cos(fc7rM)du 

1 /■* sin[(2n + l)nu] sin(/c7rw) 
sin(7ru) 



du. (A35) 



A little thought and the inequalities HA11|I and HA14|I 
show that 



lim kiTT / gn (u) sm{kinu)du = 1 — 



lim 



sin[(2ni + l)7ru] sin(fci7ru) 
sin(7rM) 



du. 



The last limit can be evaluated with the help of Eqs. (|A9I 
and l)A18|l . The result is 



lim kiTT / gn (u) sm{kiTTu)du 1 — aln 



1 + 2a 



1 - 2a 



(A36) 



The magnitude of the terms of the last sum in 
Eq. I|A33I) can be established as follows. Clearly, 



cos(fc7rt) 



E 



cos(2j7ru) 



< 



1 °° 1 

— Y - 



j=n+l 

Then, with the help of Eq. HA14|I . we have 



= 1. 



(A37) 



2sin(A:7rt) sin(2j7rt) 



■n'^ank 



j=n+l 

sin(fc7rt) 



< 



7r2a„(n + 1) 



kirt 



TTt 



< 2- 



TTt 



sin(7ri) sin(7rt) 



<7T. (A38) 



For the last two inequalities, we used the fact that 
I sin(a;)/x| < 1 as well as the fact that nt/ sin(7ri) is pos- 
itive and increasing on the interval (0, 1/2] and therefore 
it is bounded by n/2. Form Eqs. ljll3|) . (|M7jl . and (|X38|) 
we learn that 
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sup 

te(o,i/2] 



kir I gn{u) sm{k7Tu)du 





<2 + TT- 



sup 



^Oink te(0,l/2] 



sin[(2n + l)7ru] sin(fc7ru) 



< 2 + TT + 27r sup 
te(o,i/2] 



g siii(7rw) 
* sin[(2ri + l)nu] sin{kTTu) 



du 







sin(7rit) 



du 



(A39) 



where for the last inequaUty we used the relations 
7r^a„(n + 1) > 1 [see Eq. and k > n + 1. Now, 



if fc 7^ 2n + 1, Eqs. l|IT9)l and (|X39|l produce 



sup 

46(0,1/2] 



kn I gn{u) sm{kTTu)du 




< 2 + TT + (1 + tt) In 



1 + (2n+ 



1 - {2n+l)/k 

= 2 + TT + (1 + tt) In [1 + (2n + l)/fc] - (1 + tt) In [|1 - (2n + l)/fc|] 

2ra+ 1 



< 2 + TT + (1 + tt) ln(3) - (1 + tt) In 
I 



1 



k 



(A40) 



where for the last inequality we utilized the relation k > 
71 + 1 > n + 1/2 to conclude 

ln[l + (2n+ l)//c] < ln(3). (A41) 

If fc = 2n + 1, then Eqs. (|X20|l and ((X39|) produce 



sup 

te(o,i/2] 



kn j gn{u) sin{kTTu)du 





<4 + 7r + ln(2n+l). 

(A42) 



The analysis of the term 



kn / 5° (u) sin(fc7ru)du 
Jo 

is more straightforward. We start with the identity 



kiT I (7° (u) sin(fc7ru)du 
Jo 

kir f sin[(8n + 1)?™] sin(fc7ru) 
6n Jq sin(7rM) 
kn /■* sin[(2n + l)nu] sm{knu) 



sin(7ru) 

Then, from Proposition^ we readily learn that 

/■* . 1 

lim kin / g°^.{u) sm{kinu)du = 

i-^oo Jq ' ' 12a 

1 + 2a 



du (A43) 



du. 



In 



1 



1 - 8a 



In 



1- 2a 



(A44) 



Now, using the inequality fc < 4n together Eq. (|A43|I . we 



conclude 



sup 

*6(0,l/2] 



In 



;° . (u) sm{kinu)du 
1 + (8n+ l)/fc 



< 



1 + n 



1 - {8n + l)/k 



In 



1 + (2n+ l)//c 



1 - (2n+ 1)/A: 



Using Eq. (|A41|) . the fact that the function 
In [|(1 + x)/(l — a;)|] is decreasing on the interval 
(1, oo), and the fact that (8n + 1)/A: > 2, we obtain 



sup 

te(o,i/2] 



3° . (u) sin(fci7ru)du 
2 In (3) - In 



< 



1 + n 









( 


^ 2n+ 1 


)] 




k 





(A45) 



Proceeding in an analogous manner but also using 
Eq. HA20|) . one proves 



sup 

te(o,i/2] 



5° . (u) sin(A:i7ru)du 



< 



1 + n 
3 



xln(3) + - + iln(2n+l). (A46) 
3 3 

The statements of the proposition follow from 
Eqs. K36ll . (IX40ll . K42t . (IX^ . K45t . llX46ll . and the 

relation 

^n.k{t) = / [5„(m) - sin(A:7rM)dM. 
Jo 

The proof is concluded. □ 
We now focus out attention on the asymptotic behavior 
of the primitives 'ifn,k{t) on the interval (0,1/2]. We 
begin by establishing the following bound. 
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Proposition 6 Let R : {0,oo) (0,oo) be the continu- 
ous, positive, and integrable function 



R{u) 



,-1/2 



u < 1, 



w , u > 1. 

Then there is Mq > a positive constant such that 



(A47) 



r„(M) - 1 + i [gn{u) - g°{u)] 



< MoRinu) (A48) 



for all u £ (0, 1/2] and n>l. 

Proof. According to Eq. (|A3|) . we have 
7„(u, u) < u. 



(A49) 



Next, let us notice that the function sin(7ru) is positive 
and concave (convex down) on the interval (0, 1/2] and 
so, Jensen's inequality says that 

sin(7ru) = sin[(l - 2u)0 + 2M(7r/2)] 

> (1 - 2u) sin(O) + 2Msin(7r/2) = 2u (A50) 

for all u e (0, 1/2]. It follows that ii k < 2n, we have 

sin(fc7ru)^ > 

for all u < l/(4n). Using this observation, one shows 



2n 



2n 



fc— n+1 /c— n+1 

> — ^ / X ax = 



in 

for all u < l/(4n). On the other hand, if u e 
[l/(4n), 1/2], then Eqs. (|XT2|) and l|J50)l show that 

7°(u,m) = a„[l - gliu)] > a„[l - |g°(w)|] 

\ bnu J \ bn/An 

As such, for all u G (0, 1/2], the following inequality holds 



7n > ^minjl,- 



This inequality implies that 



o / \ \ 2 2 



(A51) 



for all M < 1/n. Indeed, if u < 1/n, then clearly n^u^ < 1 
and n^M^ < 56M^n^/3. Therefore, 



rt^u^ < min 1 



56u 



from which Eq. (|A51|) follows easily. Finally, from 
Eqs. and l|J5T|l . we obtain 



o / N ^ 2 2 ^ 2n nu^ nv? 



(A52) 



for all u < 1/n. Combining the last inequality with 
Eq. HA49() . one deduces 



hn{u,u) /67r2 

rn(u) = \ — < \ . 

y 7„(u, uj V nu 



Going back to Eq. ljA48|) and using the inequalities 
l5n(w)| < 1 and < 1, one obtains 



r„(u)-l + i[g„(u)-.g°(?i)] 



< k„(w)| + 2 



< \i— + = {2 + VGtt^ ] R{nu). (A53) 
" nu \ nu ^ ' 



for all u < 1/n. 

We now turn our attention to the case u> 1/n. From 
Eqs. (|Alip . (|A12|) . and (|A50(I we learn that 



\gn{u)\ < 



1 1 



< 



n sin(7rit) 2nu 



and 



1 1 



< 



3n sin(7ru) 6nw' 



(A54) 



(A55) 



respectively. Moreover, since u > 1/n, we also have 
\gn{u)\ < 1/2 and < 1/6. Let 



h{x,y) 



1-x 
1^' 



The function h{x, y) has continuous partial derivatives of 
any order on the rectangle D = [—1/2, 1/2] x [—1/6, 1/6]. 
Performing a series expansion about origin to the first 
order and employing the Taylor formula with remainder, 
we learn that for each {x,y) € D there is {x',y') g D 
such that 



h{x,y) -l + ^{x-y) = ^ 



+2 



d'^h , . .. d'^h , , ,. 



Thus, 



h{x,y) -l + -{x-y) 



< 



[x ,y) 



{x',y') 



dx^ 
dy 



A^'^y') 



dxdy 

X max{a;^,y^} < Mi max{a;^, ?/^}, 



(A56) 



where 



Ml = - sup 

^ (x',y')eD 



dx"^ 



{x ,y) 



dxdy 



{x',y') 
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is finite because the second order derivatives of h{x, y) 
are continuous on the compact set D, thus bounded. 

Combining Eq. I)A56|I with Eqs. (|A54|I and (|A55|I . we 
learn that 



r„(M) - 1 + i [gn{u) ~g°{u)] 



< Ml 



X max 



1 



1 



Ml 



for all u>l/n. The statement of the proposition follows 
immediately from Eqs. ljA53|) and (|A57ll by setting 

Afo = max |2 + V6^, Afi/4| . 

The proof is concluded. □ 

Proposition 7 Let ki oo and rii oo be two se- 
quences such that rii/ki — > a. Let 



ho{a,u) 



cos(2aui^) 



di/ 



1/2 



-1-1/2 



1 / cos(2aui/)di/ 

3 Ji 



■ cos(2aw) , 1 /"^ , 
/ di/ / cos(2aui/jdi/ 



1 

+ 2 

and define 



Ho{a) = / /io(a, u) sin(it)du. 

^0 

Lft e (0,1/2], f/ien 

lim [fc,7r«'„,,fc,(t)] = i/o(a). 

Moreover, 

sup |fc7r*„,fc(i)| < 127rAfo, 

tG(0,l/2] 

where Mq is the constant appearing in Eq. P7' 



(A58) 



(A59) 



(A60) 



Proof. We begin by performing the transformation of 
coordinates 



Jo Jo 



where LkiTTt{u) is the indicator function of the interval 
[0, kint]. As ki — > oo, we have LkiTrt{u) 1 for all t > 0. 
Let us find the limit of the sequence 



TTki 



We start with the observation that 



lim g„ 



u 

TTki 



= lim n, 

i — ^oo 



cos{2ju/ki 



J' 



= lim — 

i^oo kj 



E 

= {ni + l)/k 



E 

j=ni+l 

cos(2ui/, 



-, (A61) 



where = i/A:^. Remembering that Uijki — > a, we 
notice that the last sum is a Riemann sum approximation 
of the continuous function cosi^uv) j ^ regarded as a 
function of v over the interval (a,oo). For each z, the 
function cos(2uz/)/z/^ is evaluated on the grid made up 
of the equally spaced points {t'i,j}j>n- The mesh of the 
grid is and decays to zero a&i oo. A little thought 
shows that the limit in Eq. ljA6ip is 



lim g„ 



Trfc,: 



&v = 



cos(2auiy) 



Proceeding in a similar way, one shows 



1 

3.yi 



cos(2auz/)dj/. 



(A62) 



(A63) 



From Eqs. (|A62|I and ljA63|) the reader easily concludes 
that 



ho{a, u) sin(M) 



lim tpmM ( —77 



and therefore, 

lim IkiTTt{u)ijjni.k, 

i — ►oo 

Let us notice that 



u 
nki 



Iki7rtiu)lpni,ki ( —77 



ho (a, u) sin(M) . ( A64) 



<MoR{^]. (A65) 



The inequality is trivial if u > kiirt. Otherwise, we have 
u/{ki7T) < t < 1/2 and the inequality follows from Propo- 
sition|Sland the fact that | sin(M)| < 1. Now, the continu- 
ity of the function R(u) as well as the fact that rii/ki — > a 
imply 



limi?(^ 



= R 



Moreover, 



/■°° { niu\ ki-K f°° 

lim / R ) du = lim — ^ / R (u) du 

R [u) Au : 



TT 

a -10 



R 



au\ 

— du. 

TT / 



In these conditions, the dominated convergence theorem 
and Eq. (|A64|I allow us to conclude that 



lim 



J Iki7Tt{u)i'niM {—j^J ^"^ = J ho{a,u)sm{u)du 
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and so, the equality ljA59p is proved. 
From the inequahty HA65|I we learn that 



lim sup 

46(0,1/2] 







R (^] du 



Mo— R (u) du = Mo— < l2nMo 
n Jo n 



and the proposition is concluded. 



□ 



Proposition 8 Let ki oo and Ui — > cx) 6e two se- 
quences such that Ui/ki a and the ki 's are either odd 
or even. Let 



H{a) = Ho{a) - i<! 1 -aln 



1 + 2a 



1 - 2a 



+ 



12a 



H 


1 


+ 2a 


)-( 


1 +8a 


)] 


} 


1 


- 2a 


1 -8a 







with Ho{a) defined by Eq. iA 5. 
t e [0,1]. Then, 



Assume a 7^ 1/2 and 



lim \kiT:Gr, 



0, t = o, 

H{a), 0<t<l, 
[l-(-l)'=-]iJ(a), t=l. 

(A67) 

Moreover, there are positive constants c'l, C2, C3 and c'l 
such that 



sup \kiTGn.k{t)\ < c'l — C2 In 
te[o,i] 

if k ^ 2n + I and 



2n + 1 



k 



(A68) 



sup \k7:G„,kit)\ < 4 + In (2n + 1) , (A69) 
te[o,i] 

i/fc = 2n+ 1. 

Proof. First, let us notice that the theorem is true if 
t e [0,1/2]. Eq. (|X67jl is trivially true if t = . If t e 
(0, 1 /2], th en Eq. l|X67|) follows from Eqs. l|IT6)l . l|X30)l . 
and l|A59|l . Moreover, the existence of the constants c'/, 
C2, C3, and C4 with the desired properties follows quickly 
from Eqs. l|XT6)l . (HHJ , (HSU , and (|X60)l . and from the 
fact that G„,fc(0) = 0. 

In the remainder of the proof, we use symmetry ar- 
guments to show that the theorem holds on the whole 
interval [0, 1]. Let t > 1/2 and write 



Gn,k{t) = G'„,fc(l/2) + / [r„(u) - 1] sin(fc7ru)du 

Jl/2 

[r„(?i) — 1] sin(fc7ru)du. 

Then the symmetry relation [r„(u) — 1] sin(fc7rM) = 
(— l)'''[r„(l — u) — 1] sin[fc7r(l — u)] implies 



1 1/2 



[rn{u) - 1] sin{kTTu)du = (-l)'^G„,fe(l/2) 



and 

[r„(u) - 1] sin(fc7ru)du = (-l)'=G„,fe(l - t). 
Therefore, 

GnAi) = [1 + (-l)']G„,fc(l/2) - (-l)'=G„,fc(l - t). 
This last equation implies the inequality 

|G„,fc(t)| < 2|G„,fe(l/2)| + |G„,fc(l - t)\, 
from which it follows that 

sup |G„.fc(t)|<3 sup |G„,fe(OI- 

te(l/2,l] t6[0,l/2] 



(A70) 



>, (A66) We then have 



sup |G„,fc(t)| < max] sup |G„,fe(t)|, 
te[o,i] '-te[o,i/2] 

lup |G„,fc(t)|[<3 sup \Gn.k{t)\, 
1/2,11 J tefo,i/2l 



sup 

te(i/2,i] 



from which Eqs. l|A68p and (|A69|I follow. 

Now, \i t = 1, then G„;.fc. (1 — t) = and we have 

lim [fc,7rG„„fc,(i)] = [l + (-l)'''] lim [fc,;7rG„,.fe, (1/2)] 

= [l + (-lf']i?(a), (A71) 

where we also used the facts that the sequence {fc,;}i>o 
is either even or odd and that the theorem is true for 
< = 1/2. If 1/2 < i < 1, then < 1 - 1/2 and 

lim M [G„,.fc, (1/2) - G„,,fc, (1 - t)] =i7(a)-i/(a) =0. 

i — ^00 

Using this last observation in Eq. ljA70|l . we obtain 
lim [fc,7rG„„fc,(t)] - lim [fc,^G„,.fc, (1/2)] = i/(a). (A72) 

i^oo i—>-oo 

Eqs. (|A71|I and l|A72(l demonstrate Eq. (|A67|I and the 

proof of the proposition is concluded. □ 
We are now in the position to state the first theorem of 
Appendix A, which concerns the primitives of the func- 
tions fn,k{u). We define the functions 

Fn,k{t) = - cos(fc7ri) + fcTT [G„,fc(0 - G„,fc(l/2)] . (A73) 

Eq. (jA17p shows that the functions Fn^k/{kTr) are primi- 
tives of /„,fc. They satisfy the following theorem. 

Theorem 5 Let fc^ — > cx) and 00 be two sequences 

such that Ui/ki a and the ki 's are either odd or even. 
Assume a ^ 1/2. Then, 

lim F„,,,.,(0) = -[l + il(a)]. (A74) 

i — »-CX3 

lim F„^^kAl)^-{-l)'''[l + H{a)]. (A75) 
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If • [Oi 1] ^ M is integrable, then 



lim / i^„,,fc,(t)e(i)di = 
Jo 



(A76) 



because the integrand converges to zero pointwise by 
Eq. (|A84I) and is dominated by the integrable function 
M\^{t)\. Eq. l|J76)l foUows from Eqs. l|MT|l and (|X86)l . 

The same identities utihzed to prove Eq. ljA86|) and 
the same dominated convergence theorem also imply 



lim / F^^,k,{t) ^ - at)dt. (A77) 

IfCiu,T) : [0,1] X [0,1] — > R is integrable, then 



lim 



F„.,fc.(u)F„^,fc.(r)^'(u,T)dudr = 0. (A78) 



Finally, there are positive constants Ci, C2, C3, and C4 
such that 



2n+ 1 



sup |F„.fc(t)| < ci - C2 In 
*e[o,i] 

if k ^ 2n + I and 



sup |F„,fe(t)| < C3 + C4ln(2ra+ 1) 
te[o,i] 



(A79) 



(A80) 



ifk = 2n+l. 

Proof. Eqs. ^74)1 . (|X75|) (|X79)l, an d (|MO)l foll ow im - 
mediately from Eqs. ljX67)l . l|X68)l . llX69|) . and l|l73)l . 
One also utilizes the inequality | cos(/c7rt)| < 1. 

We turn our attention to Eq. IIA76|) . From the 
Riemann-Lebesgue lemma we learn that 



lim / ^{t) cos{k^Trt)dt = 0. 



(A81) 



For this proof, it is convenient to define the auxiliary 
function 



Glj,{t)^kTT [G„,fe(t)-G„,fe(l/2)], 
so that Eq. HA73() becomes 

Fn^kit) = - COs{kTTt) + G°„^k{t) 

Since Ui/ki a ^ 1/2, Eq. HA67P implies 
lim G°,^,^(t) = ViG(0,l). 

Also, from Eq. (|A68|1 . we learn that 
sup |G°.,,.(t)|<2c'/-24'ln^- + ^ 

te[OA] 



(A82) 
(A83) 
(A84) 

(A85) 



Since Ui/ki — > a ^ 1/2, the above equation shows that 
there are io > 1 and M > such that \G°^^j^^{t)\ < M 
for all i > io- Then the dominated convergence theorem 
implies 

hm /^(t)G:.,,,(t)dt = 0, (A86) 

»^oo Jo 



lim / ^(i)G°^,,,(t)^di = 0, 



lim / C k. (t) cos{k,7Tt)dt = 0, 



(A87) 



(A88) 



lim 



Jo 



and 



lim 



Jo 



e{u,T)Gl^^,^{u)G°„^^,^{t)dudT = 0, (A89) 



C'(u,T)G°^_fc^(u)cos(/ci7rr)dudT = 0. (A90) 



In all cases, the integrands converge to zero pointwise 
by Eq. l|X84| and are dominated by M'^\£,{t)\, M|^(i)|, 
A'P\^'{u,t)\, and M\C{u,t)\, respectively. Now, the 
Riemann-Lebesgue lemma shows that 

lim / £,{t) cosiki-Ktfdt = - / £.{t)dt 
»^°o Jo 2 7o 

+ i lim / £,{t) cos{2k,TTt)dt - [ ^{t)dt. (A91) 
The Riemann-Lebesgue lemma also implies that 



lim cos(A:i7ru) / ^'(u, r) cos(fci7rr)dr = 0. 
woo 7q 

for almost all u G [0,1] because £,'{u,t) as a function 
of r is integrable for almost all u G [0, 1] (this state- 
ment is part of the Fubini theorem). Moreover, the 
above integrand is dominated by the integrable function 
Jo ''")|d''" and then the dominated convergence the- 

orem and Fubini theorem imply 



lim 



1 fi 



Jo 



(^'(u, r) cos(/cj7ru) cos(fci7rT)dudr ~ 0. 



(A92) 

Simple algebraic manipulations by means of Eq. (jA83|) 
show that E qs. l|X87l), l|J88)l . an d (|X9T1) imply Eq. l|IT7)l , 
while Eqs. l|X89)l . l|X90)l . and llX92l) imply Eq. ^MB- 
The proof of the theorem is concluded. □ 

We conclude the section by proving the following the- 
orem, which is crucially used in the main text. 

Theorem 6 Let D C [0, 1] x [0, 1] be the triangle D = 
{{u, t) : < T < u, < u < 1}. Let K{u, t) : D ~> R be 
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a continuous function with continuous first order deriva- where the function H{a) is defined by Eq. I^A6t 
tives such that ■Ar-K(u,T) is also continuous. Then 



4n ^1 

'0 Jo 



lim n / du dTK{u,T)fn,k{u)fn,k{T) 

~ 27r2 



[1 + H{a)]^da } [K{0, 0) + K{1, 1)] (A93) 

1/4 J 

3 3 d 



[if (1, 1) - if (0, 0)] + / |-if (z., r) 



Proof. Successive integration by parts first against the 
variable r and then against the variable u produces the 
identity 



du / dTK{u,T)fn^k{u)fn^k{T) = : 



^0 



du, where 



[FnAOfK{0, 0) - 2i^„,fc(0)F„,fc(l)if (1, 0) + Fn,k{lfK{l, 1)] 



d 







d 



F„,kiur-g^K{u,u)du + Fn^kiO) I i^„,fe(u)— if(u,0)du-F„,fc(l) / F„,fc(T) — if (1, r)dT 



d 



d 



+ / F,,^k{uf^K{u,r) 







du + 


fduj 


T — U 


Jo Jo 



dTFn,kiu)Fn,kiT) K{u,t). 

OUOT 



r 



By Theorem the numbers An^k enjoy the following 
properties. Let 

M = i [|if (0, 0)1 + 2|if(l, 0)1 + |if (1,1)1] 



1 
2 



du 



|-if(u,0) 
ou 



d 



if(l,T) 











dT+ f 
Jo 




T—U 


du 



du / dr 

^0 



dudr 



K{u,t] 



Clearly, M < oo. Then from Eqs. I)A79|I and l)A80|l we 
learn that 



\An,k\ < M 



ci - C2 In 



1 - 



2n + 1 



(A94) 



ii k 2n+ 1 and 



|A„,fc| <M[c3 + C4ln(2n+l)]^ (A95) 

if A; = 2n + 1. Moreover, if fc; and are two sequences 
such that limj^oo rii/ki — a ^ 1/2 and the fc^'s are either 
odd or even, then 

lim A„^,fc^ - ^ [1 + Hia)f [if (0, 0) + if (1, 1) 

-2(-l)'=-if(l,0)] - i |-if(u,u)du (A96) 
4 Jo 



1 d 



du. 



This last property also follows from Theorem |S1 The 
straightforward argument is left for the reader. 



We now notice that the limit in Eq. (|A93() does not 
change if the summation over k is restricted to the set 
In = {k : n + 1 < k < 4:n,k ^ 2n + 1}. Indeed, the 
inequality 



\An,2n+l\ ^ nM [C3 + C4 lu (271 + 1)]' 



7r2(2n+l)2 



7r2(2n+ 1)2 



implies 



A-n.2n 



+ 1 



lim n 
n^oo 7r^(2n + l)-^ 



0, 



from which our claim follows. In these conditions, the 
limit in Eq. (|A93I) becomes 



,. ri -^-^ An k ,. ri fT-^ An.k 

hm — } — — = hm — } — — 

n — *oo TT ^ — ^ K n — *oo Tf — ^ /C 

k^Xn k^ln'^k-cvcn 



n — »oo TT * ^ iZ 



E 4^- (A97) 



We broke the sum into two parts over the even terms and 
the odd terms respectively, because these partial sums 
are easier to analyze. 

Let us begin with the sums over the even terms. We 
notice that 



■ E 

fce^Xri ;/c-Gven 



An. 



An,k _ J}_ \ ^ 

fc=[n/2] + l 



2k 
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Since 2n/(2n + 1) 1 as n — > oo, we have 



E 



1™ 2 1.2 



2n 

27r^ n^oo 2n + 1 ■'^ 



(271+1)2 



/£=[n/2] + l 

Now, for [n/2] + 1 < A: < 2n, let us define the intervals 
n + k n + k + I 



2n + l 2n + l 



By construction, the intervals are disjoint and their re- 
union is the semi-open interval 



[n/2] + 1 3n+l 



2n + l 



2n+l 



C [1/2,3/2] 



(A98) 



Now, on the interval [1/2, 3/2], we define the sequence of 
step functions 



,2k 4^2 1 



if a; e /„ 



0, otherwise. 
By the construction of the functions Cn{x), 



2n 



2n 

— y A. 

+ 1 . 



,2k- 



{2n + lf 



/2] + l 



3/2 



1/2 



In this paragraph, we use Eq. (jA96p to show that the 
sequence of functions C,i (a;) is almost surely convergent on 
the interval [1/2,3/2]. Since the function Cn{x) vanishes 
outside the interval given by Eq. (|A98|) . it follows that 
lim„^oo Cnix) = for all x < 3/4. Now, let us pick an 
arbitrary x such that 3/4 <a::<lorl<a;< 3/2. 
For each large enough n, there is a unique fc„ such that 
X e In,k„- Clearly, 



,. n + kn 
lim 

n^oo 2n -|- 1 



X, 



from which it follows that fc„/n 2x — 1. With the help 
of Eq. (|A96p . we obtain 



hm Cn {x) = lim A„ 2fe 

n^oo n — >oo 



1 



l + H 



1 



4x- 2 
K{l,l)-K{0,0) 
4(2x- 1)2 



(2n+ 1)2 

4fc2 ~ 2(2a; - 1)2 

Ki0,0) + K{1,1)~2K{1,0)] 

1 d 



2{2x - 1)2 du 



-K{u,t] 



We leave it for the reader to show that ri{x) is increas- 
ing on the interval [1/2, 1) and decreasing on the interval 
(1,3/2] (this property is not dependent on the values of 
the positive constants ci and C2). Moreover, ri{x) is in- 
tegrable on the interval [1/2, 3/2] because the singularity 
at the point x = 1 is of logarithmic type. Now, if a: G In,k 
with k > n + 1, then x > 1 and by Eq. I|A94|I we have 



\Ux)\ = 
< 477 



An,2k 

2k 



(2n-hl)2 



4fc2 



< 4|A„,2fc| 



^ ( — — 7-;— ) < 477(x). 



2n + 1 / ~ V 2n + 1 
If X e In,k with k < n, then a; < 1 and again 
2k \ . . ( n^k 



\in{x)\ <477 



2n + \ 



< At] 



2n+l 



< 4r/(a;). 



If a; is a point lying outside any /„^fc, then Cnix) = and 
the inequality Cn{x) < 4?7(x) is trivially true. 

From the facts presented in the last two paragraphs, 
we conclude that the sequence C,n{x) is dominated by 
the integrable function ■ir]{x) and is pointwise convergent 
almost everywhere to if x G [1/2, 3/4] and to 



1 



2(2a;- 1)2 



l + H 



1 



Ax -2 
1 



[if(0,0) + /^(l,l) 



.1 g 



2(2x-l)2yo du 



-K{u,t) 



du 



otherwise. An application of the dominated convergence 
theorem produces 



lim — 

n^cxi 2n 



1 

fc=[n/2] + l 
f3/2 



(2n+l)2 

i7i,2k — = hm / Cri(a;)dx 

Jl/2 



4fc2 



1 



^/3/4 2(2a;-l)2 

x[i^(o,o) + i\:(i,i)- 2/^(1,0)] 



l + H 



1 



4x- 2 

3/2 ^ 



2 ,/o 



3/4 2(2a:-^l)^ 



d 



da; 



rda; 



du 



du. 



Next, we use Eq. (|A94p to show that the sequence of 
functions C„ (a;) is dominated by some integrable function. 
Let us consider the function 



77 (a;) = M 



Performing the substitution a — 1/ (4a; — 2), we compute 

fcel„;fc-cvcn l.'^i/* ) 









Ci — C2 In ^ 


1-i 


)] 




a; 





X [K{0, 0) + K{1, 1) - 2i^(l, 0)] - ^ 

r-l g 
10 



a^(M)-mo)] + ^r|;^("'-) 



(A99) 
du. 
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The computation of the Umit 

An,k 1 



1™ ~2 ^ 12 



E 



hm 



1 



2fe-l 



2n-l 



27r^ ri^oo 2n + 1 
(2n+l)2 



fe=n+l 



(2fc- 1)2 

{2n + lf 
(2/c + l)2 



proceeds in a similar way. For [(71+ l)/2] + 1 < k < 
2n — 1, we define the intervals 



n + k n + k + I 



2n + l 2n + l 



By construction, the intervals are disjoint and their re- 
union is the semi-open interval 



,+ [(n + l)/2] + l 3n 



2n+ 1 



2ri-t- 1 



C [1/2,3/2], 



Now, on the interval [1/2, 3/2], we define the sequence of 
step functions 

{An,2k+i f2k+iyi ' if a; e In,k and k>n + l, 
An,2k~i (2fc^l)a , if a; e In,k and k < n, 
0, otherwise. 

By the construction of the functions Cn{x), 



2n+ 1 

2n-l 



E 



fc=[Il+i] + i 



(2n+l) 
{2k - 1) 



■ E ^ 

fe=n+l 



n,2fc+l 



(2n+l)2 
(2fc+ 1)2 



3/2 



1/2 



C„(a;)da;. 



I leave it for the reader to use arguments similar to 
those employed for the sums over the even coefficients 
and prove that the sequence Cnix) is dominated by the 
integrable function 4r/(a;) and is pointwise convergent al- 
most everywhere to if a; S [1/2, 3/4] and to 



1 



2(2.-1)2^^(4^) J + 

+2K{1, 0)] - ^^^-L_[x(l, 1) - if (0, 0)] 



1 d 



2(2a; 



du 



otherwise. Then, with the help of the dominated conver- 
gence theorem and of the transformation of coordinates 



a = l/(4x — 2), one concludes 



lim ^ 



fceI„;fe-odd 



B 47r2 



1/4 



x[K{fd,Q) + K{l,l) + 2K{l,Q)\ 



327r2 



X [X(l, 1) - X(0, 0)] + ^f^ ^Kiu, t) 



(AlOO) 



du. 



The statement of the theorem follows easily from 
Eqs. (jX97| . (jX99|) . and (|JTOO^ . □ 



APPENDIX B 



In this appendix, we shall prove that 



lim n^7„(u, rf = ^^^(w - r) (Bl) 



and 



lim v?%{u, t)7„(u, t) = ~ '^)' (^2) 

in the sense of distributions. More precisely, we shall 
prove the following theorem. 



Theorem 7 If h : [0, 1] x [0, 1] 

function, then 



is a continuous 



lim n 



and 



7„(m, r) r)dwdr = — j / h{u,u)du 
371"* Jo 



(B3) 



lim n' 

n— >oo 



1 /•! 

^ / 7n(w, t)7„(w, r)/l(M, T)dMdT 
Jo 



1 



Before giving the proof of the theorem, we demonstrate 
the following proposition, which will eventually allow us 
to replace 7„(m, r) by 7°(m, t) in Eqs. (|B3(I and HB4|I 
without changing the limits. 

Proposition 9 For any continuous function h{u,T), we 
have 

lim I I [7„(u,t) - 7°(u, r)]^ r)|dudr = 0. 
Jo Jo 

(B5) 

Proof. Since the function h{u,T) is continuous on the 
compact set [0,1] x [0,1], it is bounded by some posi- 
tive constant AI. Use of the inequality \h{u,T)\ < Af in 
Eq. l|B5p shows that it is enough to prove 



lim 



Jo 



[7n(u,T)-7°(M,r)]'dudr = 0. (B6) 
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To continue, we need to establish a series of useful iden- 
tities. An application of the Cauchy-Schwartz inequality 
in Eq. produces the relation 

lniu,rf < 7°(u,w)7°(r,r). (B7) 

The same Cauchy-Schwartz inequality implies 

{a + b + cf = {1-a+l-b+l-cf < 3{a^ + b^ + c^). (B8) 

On the other hand, by the orthogonality of the functions 
{sin(/c7rr), k > 1}, 

/ ln{u,TfdT V sin(A:7ru)2 ^7°(m,u). 

(B9) 



With the help of Eq. I|B8|) , one proves 

[r„(w)r„(r) - 1]' = {[r„(M) - l][r„(r) - 1] 
+ [r„(u) - 1] + [r„(r) -!]}"< 3[r„{u) - if 
X [r„(r) - if + 3[r„(ii) - 1]^ + 3[r„(r) - 1]^, 

from which it follows that 

[7„(u,t) - 7°(",t)]' = [r„(u)r„(r) - 1]^:(u,t)2 
<3[r„H~l]2[r„(r)-l]^,°(u,T)2 
+3{[r„(u)-l]2 + [r„(r)-l]2}7°(^,r)^ 



This last inequality together with Eqs. ljB7p and ljB9l 
implies 



1 /.I 



^0 



[7„(m, t) - 7° (m, r)]^ dudr 

1 . 2 

2. 0/ 



< 3 <{ y [r„(u) - 1] 7°(u,M)dM 

l2. o 



2a„ 



[r„(u) - 1] j°{u,u)du. 



Setting 



and employing the inequality q;„ < l/(7r^n) [see 
Eq. lAQIl ] , we obtain 



1 pi 



JO 



[7„(m, t) - 7°(u, r)]^ dwdr < 31^ + , I, 



This inequality shows that in order to prove the state- 
ment of Eq. HB6|I and therefore the proposition, it is 
enough to show that /„ — > as n — > oo. 



Remembering the definition of the functions r„(u), we 
compute 



Jo 

2n^/2 a„ - / \/ 7„ {u,u)y^ 7° (m, u)d-u 
"'0 

Jo 



< 



2ni/2 



1 ^ 



I- Vl-gniu)y/l'g°niu) 



where we also employed Eqs. (jA6p and l|A7(l . The func- 
tions gn{u) and satisfy the inequalities |(7n(u)| < 1 
and < 1 as well as 



\9n{u)\ < 



1 1 



n sin(7rM) 



and < 



1 1 

3?! sin(7ru) 



[for the last two inequalities, see Eqs. (|Alip and l|A12|l ] 
I leave it for the reader to prove the inequality 



for all X G [—1, 1] and a G [0, 1]. Then, clearly 



VT^VT^ > (1 - \xn (1 - \yn > 1 - Ixf - \yf 
and so 



It follows that 

2ni/2 



[\gUur + \g:{u)ndu 



< 



(1+3-'') 



1 



Q sin(7ru)" 



-du. 



Notice that the integral 



1 



-du 



Iq sin(7rM)" 

is finite for all a < 1. Choosing a = 3/4, we obtain 



/„<„-l/44fl + 3-3/4 



sin(7rM)'^/'* 



du, 



from which it follows that /„ ^ as desired. The proof 
of the proposition is concluded. □ 

Proof of the Theorem. By decomposing the function 
/i(u, r) in its positive and negative parts, which are also 
continuous functions, we may assume that h{u,T) is pos- 
itive. The proof is organized in two steps. The second 
step is left for the reader. 

Step 1. In this step, we use Proposition El to show 
that we can replace the functions 7„(u,, r) in Eqs. ljB3|l 
and (|B4|) with 7° (u, r) without changing the respective 
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limits. Because h(u, t) is positive, it turns out that it is 
convenient to consider the scalar product 

(/I.9) = / / f{u,T)g{u,T)h{u,T)dudT 
Jo Jo 

and the associated norm ||/|| = (f\f)^^'^- With the new 
notations, the claims of this step are 



lim n^||7„||' = lim n^h^ll' 



(BIO) 



and 



lim n3(7„|7„) = lim n^j°Jjr.), (Bll) 

n — ^oo n — >oo 

while the statement of Proposition El is 

lim n3||7„ - 7°||2 = or lim n^'^^W^ - = 0. 

n— >oo n— *oo 

(B12) 



Using Eq. ljB12(l and the standard inequality 

|||7n||-||7°lll<ll7n-7°ll, 

one easily argue that 

lim n3/2||^„|l = lini n^/'h°J, 



from which Eq. (|B10|) follows. On the other hand, by 
Cauchy- Schwartz inequality. 



lim 



'|(7n-7:i7n)l<(lini n3/2||^„-^°||) 



lim n^^'^W^n 



Eq. (B5) of Ref. (|T3) says that 
lim n^^'^W^n" 



(3^4)1/2 

Then Eq. I|B12|I implies 



h(u, u)du 



1/2 



< 00. 



lim n^|(7n -7^17") I = 0, 

n — ^00 

from which Eq. ljBll|l follows readily. 

Step 2. By the results of the first step, it is enough to 
prove 



lim n'^ I I 7°(u, t)^/i(m, r)dMdT = 7— 7 / h{u,u)d' 



u 

(B13) 



and 



lim n'^ 



Jo 



1n{u, t)7„(u, t)/i(u, T)dudr 



1 

— / h{u,u)du, (B14) 



47r4 

respectively. I leave it for the reader to adapt the argu- 
ments utilized in Appendix B of Ref. ifTol) and prove these 
last two statements. In the cited reference, the authors 
proved 



lim TV" ' ' 

n — *cxj 



^niu,T) h{u,T)dudT — / h{u,u)du 
Jo 37r Jq 



and the same technique can be used for Eqs. HB13|) and 
(EH- □ 
APPENDIX C: HARMONIC OSCILLATOR 

The n-th order RW-WFPI approximation of the par- 
tition function for an harmonic oscillator centered at the 
origin has the expression 

^n^iP) = / dao / dai . . . / da^nPn^iao, . . . , a4„; 
Jm Jr Jr 

where 



RW/ n\ 

p„ (ao, . . . ,a4„;/?) 



V2 



exp < —f3 



2 fl 



.k=0 



with 



1, fc = 0, 

^nAu) = { f^Afc(u), l<k<n, 
aAn^k{u), n < k < An. 

The functions Afe(M) and A„^fc(u) are defined at the begin- 
ning of Section IV. Notice that ao stands for the physical 
coordinate x. 

Let 5i^j denote the Kronecker symbol 



1, i = j, 
0, 



and define the positive definite symmetric matrix of en- 
tries 



A 



(n) _ / PrUQUJ^ 



i=J=0, 
rwise, 



\ 5ij + (3mouj^ nn,i{u)iln,j(u)du, otherwi 
with < i,j < An. Then, 

1 



P« (ao, • ■ ■ ,a4n;/3) 



:(2vr) 



and 



Zf^{l3) = 1 /[a2dct(A(")) 



1/2 



The computation of the entries of the matrix A^"^ has 
been performed numerically by Gauss-Legendre quadra- 
ture in 8n points. 
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